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1.1. Wstęp — częstościowa interpretacja prawdopodobieństwa 

W otaczającej nas rzeczywistości mamy do czynienia nie tylko z doświadczeniami 
deterministycznymi ale też z tak zwanymi doświadczeniami losowymi, których wyniku nie 
potrafimy przewidzieć, jednak potrafimy określić zbiór możliwych wyników. Rachunek 
prawdopodobieństwa zajmuje się opisem doświadczeń losowych, które można powtarzać w 
tych samych warunkach oraz badaniem prawidłowości którym podlegają. Zaczniemy od 
opisu pewnego doświadczenia wykonanego przez Rutherforda i Geigera. Obserwowali oni 
emisję cząstek radioaktywnej substancji w przedziałach czasu o długości 7,5 sek. Dokonano 
n = 2608 doświadczeń. Zaobserwowane wyniki przedstawia tabela, gdzie k oznacza liczbę 
cząstek wyemitowanych w przedziale czasu o długości 7,5 s, n(k) oznacza liczbę obserwacji 


10 
w których wystąpił wynik k, 2 n(k) =n, w(k) jest częstością pojawienia się wyniku k, 


k=0 
w(k)= mA ; 
n 
k |n(k) w(k) 
0 33 0,022 
1 203 0,078 
2 383 0,147 
3 525 0,201 
4 532 0,204 
5 408 0,156 
6 273 0,105 
7 139 0,053 
8 45 0,017 
9 27 0,010 
210 |16 0,006 


Wartości funkcji w(k) dla innej serii pomiarów mogą być inne. 


Zauważono, że dla długich serii pomiarów wartości funkcji w(k) stabilizują się wokół 
—A 9k 


wartości funkcji p(k) = > gdzie A jest pewną stałą (w podanym przykładzie 


k! 
A=3,877871). Wartości funkcji p(k) nie zależą od n. 


k_ |n(k) w(k) p(k) 
0 53 0,022 0,021 
1 203 0,078 0,081 
2 383 0,147 0,156 
3 525 0,201 0,201 
4 532 0,204 0,195 
5 408 0,156 0,151 
6 273 0,105 0,097 
7 139 0,053 0,054 
8 45 0,017 0,026 
9 27 0,010 0,011 
210 |16 0,006 0,007 
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liczba emitowanych cząstek 


Inne przykłady doświadczeń losowych to rzut kostką, rzut monetą, obserwacja czasu 
niezawodnej pracy pewnego urządzenia,.. 


1.2. Definicja i własności prawdopodobieństwa. 

Zbiór wszystkich możliwych wyników doświadczenia losowego nazywamy przestrzenią 
zdarzeń elementarnych i oznaczamy przez Q .Elementy zbioru Q nazywamy zdarzeniami 
elementarnymi i oznaczamy przez 0. 


Przykłady przestrzeni zdarzeń elementarnych. 
1. W doświadczeniu Rutherforda i Geigera obserwowano liczbę cząstek wyemitowanych 
przez ciało radioaktywne - Q = {0,1,2,....}. 
. Wykonujemy rzut kostką - obserwujemy liczbę oczek na kostce - Q =(1,2,3,4,5,6) 
Obserwujemy czas £ niezawodnej pracy urządzenia - (2=ft:te [0,oo)|. 
4. Wyznaczamy czas ściągania pliku i liczbę przerw w jego ściąganiu - 
0=[(t,k):te [0,»),ke {0,1,2,....}}. 
5. Rzucamy symetryczną monetą aż do pojawienia się orła - Q =[O,RO,RRO,....). 


© N 


Zdarzeniem losowym nazywamy podzbiór Q. Zbiór wszystkich zdarzeń oznaczać będziemy 
literą F. Dla zdarzeń, podobnie jak dla zbiorów, określamy alternatywę (sumę) zdarzeń, 
koniunkcję (iloczyn) zdarzeń, różnicę zdarzeń. Zdarzeniem przeciwnym do zdarzenia A 
nazywamy zdarzenie A =Q\A. Q nazywamy zdarzeniem pewnym, Ø - zdarzeniem 
niemożliwym. Jeżeli Af B = 9, to mówimy, że zdarzenia A, B wykluczają się. 

Nie wszystkie podzbiory Q będziemy uważać za zdarzenie i zaliczyć do zbioru F. 


Ze względu na wymagania dotyczące działań na zdarzeniach, zakładamy, że zbiór zdarzeń F 
spełnia następujące warunki: 
L.F źó, 


2A€F>AEF, 


3.4,,4,,...€ F >|JA,eF 
n=l 
Rodzinę F nazywamy 6-ciałem zbiorów. 
Z podanych warunków wynika, że Q=AvAeF,jeśli Ae F oraz, $=QeF. 


W przypadku dyskretnej (tzn. skończonej lub przeliczalnej) przestrzeni Q zdarzeniem 
losowym może być dowolny jej podzbiór. W przypadku nieprzeliczalnej przestrzeni Q 
rodzinę F należy precyzyjnie określić. W zagadnieniach dobrym modelem probabilistycznym 
okazują się pewne podzbiory prostej, płaszczyzny lub przestrzeni. Za rodzinę F przyjmuje się 
wówczas o-ciało zbiorów borelowskich na prostej, płaszczyźnie lub w przestrzeni, przez co 
rozumie się najmniejsze ©-ciało zbiorów zawierające przedziały otwarte (a, b) w R (koła 
otwarte w R*, kule otwarte w R”) Mówiąc obrazowo, zbiór borelowski w R (w R* „Kk=2,3) to 
każdy zbiór, który można otrzymać jako wynik przeliczalnych działań mnogościowych 
wykonanych na rodzinie wszystkich przedziałów na prostej czy kul otwartych w przestrzeni 
k-wymiarowej (k=2,3). Na przykład zbiorem borelowskim w R jest każdy przedział 
jednostronnie czy dwustronnie domknięty i każda półprosta. 


Opisując doświadczenia losowe podajemy zbiór wszystkich możliwych wyników, określamy 
rodzinę F ale najbardziej interesujące są pytania o szansę zajścia poszczególnych zdarzeń. Tę 
szansę określającą zajście zdarzenia mierzy funkcja zwana prawdopodobieństwem. 


Definicja 1.2. 
Prawdopodobieństwem nazywamy funkcję określoną na rodzinie zdarzeń F spełniającą 
następujące warunki: 

1.0<P(A)<1, dla każdego AeF 

2.P(Q)=1 

3. Jeżeli A; ,A,,...są parami rozłączne (tzn. A, © A j= é, dla dowolnych i+ j), to 


PIJA) =Ñ PA,) . 


n=l n=l 


W szczególności P(AUB)= P(A)+ P(B) dla zdarzeń rozłącznych A, B. 


Z powyższych warunków (aksjomatów) wynika wiele własności prawdopodobieństwa, 
z których najważniejsze są przedstawione poniżej: 


Fakt 1.3. 
1. Jeżeli ACB, to P(A)< P(B), 
2. P(A?) = 1- P(A), a stąd P($) =0. 
3. Jeżeli ACB, to P(BNA)=P(B)-P(A)=P(B)-P(A0B) , 
4. P(AUB)=P(A)+P(B)- P(ANB). 


5, PIJA) < Y P(A,) 
i=l i=1 


6. Jeżeli A, C A, ©..,to P((_]JA,)= lim P(A,), 


n=l 


7. Jeżeli A, > A, >...to P(f)A,)=limP(A,). 


n=l 


W przypadku trzech zbiorów A, B, C wzór z punktu 4. przyjmie postać: 
P(AUBUC)=P(A)+P(B)+P(C)-P(ANB)-P(BNC)-P(CNA)+P(ANBNC) a 
przypadku dowolnej, skończonej liczby zbiorów wzór ten można uogólnić w następujący 
sposób: 


P(JA, = È P(A)- ),P(Ą, NA) +...+(-1)"* P(A 0...0A,) 


i=l 1<iśn Ii, <i, Sn 


Fakt.1.4. 
Jeżeli Q = {ø@,; :ie I}, gdzie I jest zbiorem skończonym lub /=N oraz P(0,)=p,, przy 


czym p, 20 oraz zp =1, todla A c Q wzór P(A)= Èp, określa 


i=1 0;,EA 


prawdopodobieństwo na rodzinie wszystkich podzbiorów zbioru Q. 


Trójkę (Q, F, P) nazywamy przestrzenią probabilistyczną. 
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e Przykład 1.1. 
Niech A, B,C będą zdarzeniami. Zapisać za pomocą działań na zbiorach następujące 
zdarzenia: 
a) zachodzi dokładnie jedno ze zdarzeń A, B, C; 
b) zachodzą dokładnie dwa spośród zdarzeń A, B,C; 
c) zachodzą przynajmniej dwa spośród zdarzeń A, B,C. 
d) zachodzą co najwyżej dwa spośród zdarzeń A, B,C. 
Rozwiązanie. 
a) Zachodzi dokładnie jedno ze zdarzeń A, B, C, to oznacza dokładnie, że zachodzi A i nie 
zachodzą B ani C lub zachodzi B i nie zachodzą A ani C lub zachodzi C i nie zachodzą 
A ani B czyli, że zachodzi zdarzenie 


AME NC UG NBENCUGNETTE. 


b) Podobnie — zachodzą dokładnie dwa spośród zdarzeń A, B,C oznacza, że zachodzi 
zdarzenie 


ANBNCUANBNCUANBNC. 


c) Zachodzą przynajmniej dwa spośród zdarzeń A, B,C oznacza, że zachodzą dokład- 
nie dwa spośród zdarzeń A, B,C lub zachodzą wszystkie trzy zdarzenia, czyli zachodzi 
zdarzenie 


ANBNCUANBNCUANBNCUANBNC. 


d) Zachodzą co najwyżej dwa spośród zdarzeń A, B,C oznacza, że nie zachodzą wszystkie 
trzy zdarzenia, czyli zachodzi zdarzenie 


(ANBOREV=4UB' UGC 


e Przykład 1.2. 
Studenci Wydziału Elektroniki muszą zaliczyć dwa lektoraty: z języka angielskiego i z ję- 
zyka niemieckiego. Z danych Dziekanatu wynika, że 3 studentów zalicza lektorat z języka 
angielskiego, oba lektoraty zalicza co czwarty student, zaś przynajmniej jeden z lekto- 
ratów zalicza również Z studentów. Jakie jest prawdopodobieństwo, że losowo wybrany 
student: 

a) nie zaliczył żadnego lektoratu? 

b) zaliczył język angielski i nie zaliczył języka niemieckiego? 

Rozwiązanie. 

Niech A oznacza zdarzenie ”losowo wybrany student zaliczył lektorat z języka angielskie- 

go”, przyjmujemy, że P(A) = 2, B - zdarzenie ” losowo wybrany student zaliczył lektorat 

z języka niemieckiego”. 

a) Oczywiście chodzi o zdarzenie Æ’ N B’, więc 
P(ANnNB')=P((AUB)/)=1-P(AUB)=1-— 3 = 


1 
3 


b) Podobnie: P(ANB) = 4, P(AUB) = 2, więc P(A\B) = P(4A)-P(ANB) = 2— 


1 

4 

e Przykład 1.3. 
Studenci Wydziału PPT zdają w sesji zimowej I roku egzaminy z przedmiotów A,B,C. 
Wiadomo, z danych poprzednich lat, że przedmiot A zalicza 60% studentów, przedmiot B 
zalicza 80% studentów i przedmiot C zalicza 10% studentów. Studenci, którzy zaliczyli A 
i B stanowią 55% ogółu, ci którzy zaliczyli A i Č stanowią 45% ogółu a studenci, którzy 
zaliczyli B i C stanowią 60% ogółu. Sesję zimową zalicza ok. 40% studentów. Obliczyć 
prawdopodobieństwo, że losowo wybrany student: 


a) zaliczył przynajmniej jeden egzamin, 
b) zaliczył przynajmniej dwa egzaminy. 
Rozwiązanie. 


Wiemy, że P(A) = ł, P(B) = %4, P(C) = 5 oraz 


P(ANB) =, P(ANC) =, P(CNB)=<5, P(ANBNC) = 


100? 100? 100 100: 


Zatem: 

a) P(AUBUC) = P(A)+P(B)+P(C)-P(ANB)-P(BNC)-P(CNA)+P(ANBNC) = 
tito i e 16 P 6 = 15 

b) P(ANBUANCUBNO)) =P(AnNB)+P(ANC)+P(BNC) - P((ANB)nN(AN 
C)) - P((ANB)N(BNC)) - P((ANC)nN(BNO))+P((ANB)N(ANC)N(BNC)) = 
P(ANB) + P(ANC) + P(BNC)-2P(ANBNC)=4,. 


1.1 Przykłady przestrzeni probabilistycznych. 


1.1.1 Prawdopodobieństwo klasyczne 


Przestrzeń Q = fwj,wa,...,w„) jest zbiorem n zdarzeń elementarnych, z których każde 
zachodzi z tym samym prawdopodobieństwem, czyli P(łwy)) = + dla k = 1,2,... n. 
Zgodnie z Faktem 1.4. wzór 


i (1) 


gdzie A oznacza liczbę elementów zbioru A, określa prawdopodobieństwo na wszystkich 


zdarzeniach A C Q. Jest to tzw. prawdopodobieństwo klasyczne. 


W rozwiązywaniu zagadnień, w których przestrzeń zdarzeń elementarnych jest skończona 
przydadzą się nam wiadomości z kombinatoryki. 


Podstawowe schematy kombinatoryczne 
Niech A oznacza dowolny zbiór n różnych elementów A = f0,,0,...,an). 
e Wariacje z powtórzeniami. 


k-wyrazową wariacją z powtórzeniami zbioru A nazywamy każdy k-wyrazowy ciąg 
elementów tego zbioru Liczba V;Ę wszystkich wariacji k-wyrazowych z powtórzeniami ze 
zbioru n-elementowego wynosi 


Wam (2) 


e Wariacje bez powtórzeń. 


k-wyrazową wariacją bez powtórzeń zbioru A (k < n)nazywamy każdy k-wyrazowy 
ciąg różnych elementów tego zbioru Liczba V; wszystkich k-wyrazowych wariacji bez 
powtórzeń ze zbioru n-elementowego wynosi 


VE = n(n = 1)(n-2)::::- (n=k+1) (3) 


Jeżeli k = n, to k-wyrazową wariację bez powtórzeń ze zbioru A nazywamy n-wyrazową 
permutacją. Zatem liczba P„ wszystkich permutacji zbioru n-elementowego wynosi 


Pm nh: (4) 
e Kombinacje. 


k-elementową kombinacją bez powtórzeń z n-elementowego zbioru A nazywamy 
każdy k-elementowy podzbiór zbioru A. Liczba C$ wszystkich kombinacji k-elementowych 
bez powtórzeń ze zbioru n-elementowego wynosi 


C: = (i). (5) 


e Kombinacje z powtórzeniami. 


k-elementową kombinacją z powtórzeniami z n-elementowego zbioru A nazy- 
wamy zbiór składający się z k elementów, różnych lub nie różniących się między sobą, 
wybranych spośród n elementów, przy czym nie jest istotna kolejność, w jakiej elementy 
tego zbioru są ustawione. Liczba G wszystkich k-elementowych kombinacji z powtórze- 
niami ze zbioru n-elementowego jest taka sama, jak liczba wszystkich różnych zbiorów 
k-elementowych, jakie można utworzyć ze zbioru (n+-k-1) elementów. Kombinacje z po- 
wtórzeniami można interpretować, jako rozmieszczenie k nierozróżnialnych elementów w n 
szufladach, przy czym w każdej szufladzie może być dowolna, nie przekraczająca k, liczba 


elementów. 
= k+n-1 k+n=i 
R Ja) © 


e Przykład 1.4. 


W teorii cząstek elementarnych bada się rozmieszczenie k cząstek w podzielonej na komór- 
ki przestrzeni fazowej, którą można matematycznie opisać np.jako podzbiór przestrzeni 
czterowymiarowej, gdzie współrzędnymi są położenie i pęd cząstki. Fizycy stwierdzili do- 
świadczalnie, że niektóre cząstki zachowują się, jak kule rozróżnialne, inne - jak kule 
nierozróżnialne i zaproponowali trzy następujące modele zachowania się cząstek : 

a) statystyka Maxwella-Boltzmanna. Cząstki zachowują się, jak kule rozróżnialne, więc 
pytając o liczbę możliwych rozmieszczeń k cząstek w n komórkach mamy do czynienia z 
wariacjami z powtórzeniami i każde spośród n* rozmieszczeń jest jednakowo prawdopo- 
dobne. Nie znaleziono jeszcze cząstek, które zachowywałyby się zgodnie z tym modelem. 


b) statystyka Fermiego-Diraca. Cząstki zachowują się, jak kule nierozróżnialne, ale w każ- 
dej komórce może być co najwyżej jedna cząstka i wszystkie możliwe rozmieszczenia są 
jednakowo prawdopodobne. Tak zachowują się np. elektrony, protony i neutrony. 


c) statystyka Bosego-Finsteina. Cząstki zachowują się, jak kule nierozróżnialne, w każdej 


komórce może być dowolna liczba cząstek i wszystkie możliwe rozmieszczenia są jednako- 
wo prawdopodobne. Tak zachowują się np. fotony. 


Zadanie — w każdym z rozważanych wyżej modeli wyznaczyć prawdopodobieństwo, z 
jakim k (k < n) cząstek można rozmieścić po jednej w k z góry zadanych komórkach. 
Rozwiązanie. 
a) Przy ustalonej permutacji k rozróżnialnych cząstek tylko jedno rozmieszczenie spośród 
wszystkich ně możliwych rozmieszczeń spełnia żądany warunek. Ponieważ cząstki moż- 
na ustawić na k! sposobów, więc prawdopodobieństwo, z jakim k rozróżnialnych cząstek 
można rozmieścić po jednej w k z góry zadanych komórkach równe jest E. 
b) Jeżeli nie odróżniamy cząstek, ale w każdej komórce może być co najwyżej jedna to 
wystarczy wybrać k spośród n komórek i wrzucić do niej cząstkę, a to można zrobić na 
,) sposobów. Tylko jeden z nich spełnia warunek z zadania, więc prawdopodobieństwo, 
z jakim k nierozróżnialnych cząstek można rozmieścić po jednej w k z góry zadanych 


komórkach wynosi T 
k 


c) W tym przypadku mamy do czynienia z k-elementowymi kombinacjami z powtórze- 
niami z n różnych elementów. Dla ustalonych k-komórek jest tylko jeden ciąg cząstek 
spełniający warunki zadania, więc prawdopodobieństwo, z jakim k nierozróżnialnych czą- 
stek można rozmieścić po jednej w k z góry zadanych komórkach wynosi TT 


e Przykład 1.5. 


W pudle są kule białe i czarne. Razem jest ich n. Ile powinno być kul czarnych, aby praw- 
dopodobieństwo wylosowania (bez zwracania) dwu kul różnych kolorów było takie samo, 
jak prawdopodobieństwo wylosowania dwu kul tego samego koloru? 


Rozwiązanie. 


n 


Dwie spośród n kul można wybrać na ( ) sposobów. Oznaczmy przez k liczbę kul czar- 


2 
nych. Zdarzeniu A ”wylosowano dwie kule różnych kolorów” sprzyja k(n — k) zdarzeń 


elementarnych. Zdarzeniu B ” wylosowano kule tego samego koloru” sprzyja (o) + S 
zdarzeń elementarnych. Wykorzystując wzór na prawdopodobieństwo klasyczne otrzymu- 
jemy 

P(A) — k(n=k)-2 Śr P(B) e 2k—2nkpn>n. 


n(n—1) n(n—1) 


Ponieważ zdarzenia A i B są przeciwne, to zamiast warunku P(A) = P(B) wystarczy 
rozważać jeden z warunków P(A) = 3 lub P(B) = 3. Każdy z nich jest równoważny 
równaniu 

4k? — 4kn +n? —n=0. 
Rozwiązaniami tego równania są liczby kų = YB. oraz k= mA, Zauważmy, że jeżeli 
yn nie jest liczbą naturalną, to zadanie nie ma rozwiązania. Jeżeli zaś y/n jest liczbą 


naturalną, to zarówno 


— ny _ yava) — ntyń _ vaat) 
M ERIE WEGO O REI WERDER 


jak i 
są liczbami naturalnymi oraz kı + ko = n. 


Podsumowując — zadanie ma rozwiązanie jedynie w przypadku, gdy n jest kwadratem 


liczby naturalnej (tylko wówczas „yn jest liczbą naturalną), czarnych kul powinno być 
k= AL — VaI) Jub ko = ntyn — vevan) 


2 2 g 


e Przykład 1.6. 
W szufladzie są dwie skarpety na prawą nogę i jedna na lewą nogę. Prawdopodobieństwo, 
że losowo wybierając dwie skarpety otrzymamy parę równe jest 5 = 2, zaś prawdopodo- 
2 


() 


bieństwo wyciągnięcia dwu prawych wynosi (e) = 3. Do szuflady dołożono jedną skarpetę. 
2 


Jaka to jest skarpeta, skoro teraz prawdopodobieństwo, że wylosowane dwie skarpety sta- 
nowią parę, wynosi >? 

Rozwiązanie. 

Wykorzystajmy poprzedni przykład. Mamy n = 4. Wylosowanie pary skarpet odpowia- 
da wylosowaniu kul różnych kolorów. Zatem skarpet jednego typu może być kı = 3 lub 
ką = 1, czyli dołożono prawą skarpetę. 


e Przykład 1.7. 


Przypuśćmy, że do jeziora zawierającego nieznaną liczbę N ryb wpuszczono dodatkowo 
1000 ryb oznakowanych (np. pomalowanych na czerwono). Po pewnym czasie dokonano 
połowu 1000 ryb i znaleziono wśród nich 100 ryb z czerwonymi plamami. Jak na podsta- 
wie tych danych ocenić liczbę ryb w jeziorze? 

Rozwiązanie. 

Za ocenę N przyjmiemy taką liczbę, dla której prawdopodobieństwo wyłowienia 100 zna- 
czonych ryb spośród 1000 jest największe. Przyjmując, że liczba ryb w jeziorze równa 
się (N + 1000), wyznaczymy prawdopodobieństwo py(4), gdzie A oznacza zdarzenie po- 
legające na wylosowaniu 100 ryb oznaczonych przy losowaniu 1000 ryb. Q jest zbiorem 
kombinacji 1000-elementowych ze zbioru N + 1000 elementowego. Wśród wylosowanych 
jest 900 nieoznakowanych, więc N > 900. Stąd 


©G_ (N+1000 = _ (1000)/ N . — (o) ooo) 
= ( 1000 ), A= (io) (l: więc pn(4) = (BED) . 
Aby określić najbardziej prawdopodobną liczbę ryb w jeziorze, wyznaczymy wartość N, 
przy której py(A) osiąga wartość maksymalną. Rozpatrzmy 
PN(A)_ _ N2 — | + z 100N+900-1000 
PN-1(4) _ (N-900)(N+1000) (N —900)(N +1000) ` 
Zauważmy, że iloraz ten jest większy od 1 bądź mniejszy niż 1 w zależności od tego, czy 
100N < 900 - 1000, czy 100V > 900 - 1000. Oznacza to, że gdy N rośnie, liczby py(4) 
najpierw rosną a potem maleją. Rozważany iloraz osiąga wartość największą, gdy N jest 


największą liczbą naturalną nie przekraczjącą 200-1000, czyli N = 2001000 lub N = 8999. 


1.1.2 Przeliczalna nieskończona przestrzeń probabilistyczna 


Niech Q = {w1,w2,... }. Jeżeli P(w;) = pi przy czym p; > O oraz 5) p; = 1, to prawdopo- 
i=1 
dobieństwo zdarzenia A C Q określone jest wzorem 


P(AJ= > Pi. (7) 


e Przykład 1.8. 
Dwaj gracze, A oraz B, rzucają na przemian monetą, dopóki dwa razy pod rząd upadnie 
ona na tę samą stronę. Jeżeli drugi pod rząd orzeł albo druga pod rząd reszka pojawi się 
w rzucie nieparzystym, to wygrywa gracz A. W przeciwnym przypadku wygrywa gracz 
B. Obliczyć prawdopodobieństwo wygranej dla każdego z graczy. 
Rozwiązanie. 
W opisanej grze zdarzeniem elementarnym jest ciąg, którego elementami są orły lub reszki 
i na ostatnich dwu miejscach, po raz pierwszy pod rząd są dwa orły lub dwie reszki, czyli 
Q = foo,rr,orr,roo,oroo,rorr,...| ma nieskończenie, ale przeliczalnie, wiele zdarzeń 
elementarnych. Niech a; oznacza zdarzenie ” druga reszka pod rząd pojawiła się po raz 
pierwszy w k-tym rzucie monetą”, zaś 6% - ”drugi orzeł pod rząd pojawił się po raz 
pierwszy w k-tym rzucie.”. Oczywiście k = 2,3,.... Na przykład as = ororr, (5 = roroo. 
Prawdopodobieństwo dla zdarzeń ay, 8; wynosi 


P({ar}) = PR) = 2 k=2,3,.... 


Gracz A wygra, jeżeli zajdzie zdarzenie A = faz, 63, a5, 85,...|. Wygranej gracz a B 
sprzyja zdarzenie A = 4ao, 62, 04, 64,...|. Zatem 


oraz 


P(B) = È P({ox fa} = 227% = 2-22. aR. 


Prawdopodobieństwo wygranej gracza, który rzuca monetą na parzystych miejscach jest 
dwa razy większe niż gracza, który rzuca monetą na nieparzystych miejscach. W tej grze 
pozwólmy przeciwnikowi rozpocząć grę! My rzucajmy na miejscach parzystych! 


1.1.3 Prawdopodobieństwo geometryczne 


W wielu zagadnieniach zbiór Q można przedstawić jako podzbiór borelowski przestrzeni 
R (R? lub R”), który ma skończoną miarę m - długość (pole, objętość). Dla zdarzenia 
losowego A C Q, dla którego tę miarę potrafimy obliczyć, prawdopodobieństwo zdarzenia 
A zdefiniowane wzorem 


P(A) = ——, 8 
nazywamy prawdopodobieństwem geometrycznym. 


e Przykład 1.9. 
Kawałek drutu o długości 20cm zgięto pod kątem prostym w przypadkowo wybranym 
punkcie. Następnie zgięto drut jeszcze w dwu punktach tak, by powstała ramka prosto- 
katna o obwodzie 20cm. 
a) Jakie jest prawdopodobieństwo, że pole ograniczone ramką nie przekroczy 21cm? ? 
b) Jakie jest prawdopodobieństwo, że pole ograniczone ramką jest równe 21cm?? 
Rozwiązanie. 
a) Niech z oznacza odległość wybranego punktu od bliższego końca drutu. Wówczas 
Q = [0,10]. Zdarzenie A ” pole ograniczone ramką nie przekracza 21cm?” zachodzi wtedy 
i tylko wtedy, gdy 1(10 — x) < 21. Rozwiązując nierówność 


—g” +10x — 21<0 dla z€([0,10] 


otrzymujemy 


A = {x € [0,10] : z € [0,3] U [7,10]}, więc P(A)= 28 =$. 


Z 


b) Niech B oznacza zdarzenie ”pole ograniczone ramką jest równe 21cm”. Wówczas B 
zachodzi wtedy i tylko wtedy, gdy z = 3 lub z = 7, więc P(B) = ae = 0, ponieważ 
długość zbioru złożonego z dwu punktów wynosi 0. Zauważmy, że zdarzenie B jest możliwe 


ale prawdopodobieństwo jego zajścia równe jest 0. 


e Przykład 1.10. 


W każdej chwili odcinka czasu T jednakowo możliwe jest nadejście do odbiornika każde- 
go z dwu sygnałów, które w tym odcinku czasu zostaną przesłane. Odbiornik nie może 
przyjąć drugiego sygnału, jeżeli nadejdzie on w czasie krótszym niż T od chwili nadejścia 
pierwszego sygnału. Obliczyć prawdopodobieństwo przyjęcia przez odbiornik obu sygna- 
łów. 

Rozwiązanie. 

Niech x i y oznaczają czasy nadej- 


A 
ścia sygnałów do odbiornika. Wtedy Y 
przestrzeń zdarzeń elementarnych T 
Q= {(x,y): x,y € [0,T]} 
możemy interpretować jako kwadrat P 
o boku T a interesujące nas zdarze- 
nie można zapisać w postaci O > T fi 


A={(xz,yY ETxT:|r-ylļ| 2T}. 


co 


Rys. 1.01. 
Zatem P(A) 


e Przykład 1.11. 


Z przedziału [0,1] wybieramy losowo trzy liczby x,y,z. Jakie jest prawdopodobieństwo, 
że ich suma jest liczbą z przedziału |5, 1]? 

Rozwiązanie. 

W tym przykładzie Q = { (x,y,z): 0 < x,y,z < 1}, czyli geometrycznie Q jest sześcianem 


jednostkowym. Rozważane zdarzenie to zbiór 
A= Hez) : 4<r+y+z< i: 


Tutaj m(4) jest objętością zbioru A, który jest różnicą dwu ostrosłupów. Zatem 
m(A) = 3(31:1-4-4:5)=4:3=4. 


Ponieważ m(9) = 1, więc P(A)= £. 


2 Prawdopodobieństwo warunkowe. 
Niezależność zdarzeń. 


2.1 Prawdopodobieństwo warunkowe. 


Niech (Q, F, P) będzie przestrzenią probabilistyczną, a B - dowolnie ustalonym zdarze- 
niem takim, że P(B) > 0. 


Definicja 2.1. Prawdopodobieństwem warunkowym zdarzenia A pod warunkiem B nazy- 


wamy liczbę 


PLA|B) = Ery” 


Stąd oczywiście P(AN B) = P(A|B)P(B). 


Posługując się zasadą indukcji matematycznej możemy udowodnić, że dla dowolnego n, 
przy założeniu, że P(A N --- N An-1) > 0, prawdziwa jest równość 


PAren A = PANPA NAA 042) +. PA. ANNA). (9) 


Fakt 2.2. Jeżeli P(B) > 0, to funkcja P(.|B) określona na F spełnia aksjomaty praw- 
dopodobieństwa. 


Przykład 2.1. 
Rzucamy dwa razy symetryczną kostką. 
a) Jakie jest prawdopodobieństwo wyrzucenia różnej liczby oczek? 


b) Jakie jest prawdopodobieństwo wyrzucenia różnej liczby oczek, jeżeli suma oczek wy- 
nosi 11? 

c) Jakie jest prawdopodobieństwo wyrzucenia różnej liczby oczek, jeżeli suma oczek wy- 
nosi 10? 

Rozwiązanie. 

Przestrzeń zdarzeń elementarnych Q jest zbiorem par uporządkowanych (a, b), gdzie 
a,be41,2,3,4,5,6). Wszystkie zdarzenia elementarne są jednakowo prawdopodobne. 


a) Niech A oznacza zdarzenie ” wypadła różna liczba oczek”, czyli 
A= {(a,b) EQ: a+b} = Q \ {(1, 1), (2,2), (3,3), (4,4), (5,5), (6,6)4. 


Ponieważ 0 = 6? = 36 oraz A = 36 — 6 = 30, więc P(A) = ¥ = £. b) Zdarzenie ” suma 
oczek wynosi 11” oznaczmy przez B. Oczywiście B = £(5,6),(6,5)|. Ponieważ B C A, 
więc BN A= B i stąd 


— 


W tym przykładzie informacja o zdarzeniu B dawała pewność, że zajdzie zdarzenie A. 
c) Warunkiem jest zdarzenie D = +(4,6), (6,4), (5,5)). Mamy: 


P(D)=+2, AND ={(4,6),(6,4)}, P(AND)=L. 


Zatem 


P(A|D) = "p = 3. 


Jak widać, przy różnych warunkach prawdopodobieństwo zajścia zdarzenia A jest różne. 


© Przykład 2.2. 
Wybrano losowo dwie liczby z przedziału [0,1]. Jakie jest prawdopodobieństwo, że zy > 
0,09 , jeżeli wiadomo, że z+y<1? 
Rozwiązanie. 


Przestrzeń zdarzeń elementarnych jest kwadratem jednost- 4” 
kówyń Q = {(x,y): x,y € [0,1]}. 1h 
Interesujące nas zdarzenia to A= {(x,y): zy > 0.09) 
oraz B= (ay): x+y < 1}. Oczywiście P(B) = 4. 
Ponieważ 


ANB = {(z,y): 01<v<0,9; 2 <y<1- z}, 


więc P(AN B) = f Q- x- =P)da = 2 — 7% ln9. 
0,1 


Stąd P(A|B) = FG) 


I 
NO 
TT" 
GN 
| 

Ę 
— 
D 
© 

— 


e Przykład 2.3. 


Studenci Wydziału Elektroniki muszą zdać w I semestrze trzy egzaminy: z fizyki (A), 
analizy matematycznej (C) i z algebry (B). Z danych Dziekanatu wynika, że 10% stu- 
dentów zalicza I semestr a 90% - zdaje egzamin z fizyki. Jeżeli student zaliczy algebrę i 
fizykę, to prawdopodobieństwo, że zda analizę wynosi É, Jakie jest prawdopodobieństwo, 
że student, który zdał fizykę, zda algebrę? 

Rozwiązanie. 


Skorzystamy ze wzoru 
P(AN BAC) = P(A); P(B|A)- P(C|AN B). 
Mamy 5 = 4: P(B|A): $, skąd P(B|A) = $. 


I0 10` 
2.2 Wzór na prawdopodobieństwo całkowite i wzór Bayesa. 


Załóżmy,że Q jest sumą rozłącznych zbiorów B; € F dla i€. Wówczas dla dowolnego 


zdarzenia A zbiory A N B; są parami rozłączne. Ponadto A= U (AN B;), więc P(A)= 
iEI 

> P(ANB;). Jeżeli wszystkie zdarzenia B; mają dodatnie prawdopodobieństwo, to P(AN 

iE 


B,)=P(A|B;) - P(B;) dla każdego ie] i otrzymujemy następujące twierdzenie. 


e Twierdzenie 2.1. (Twierdzenie o prawdopodobieństwie całkowitym) 


Jeżeli Q jest sumą rozłącznych zbiorów B;, przy czym P(B;,) > 0 dla wszystkich tel, to 
dla dowolnego zdarzenia A zachodzi równość 


P(A)= > ,P(A|B;) - P(B,). 


icl 


Czasem zdarzenia B; występujące we wzorze na prawdopodobieństwo całkowite nazywa- 
my przyczynami, zdarzenie A - skutkiem. Rozważmy zagadnienie w pewnym sensie 
odwrotne do zagadnienia obliczania prawdopodobieństwa całkowitego. Mianowicie zapy- 
tajmy, jakie jest prawdopodobieństwo przyczyny B;, gdy znany jest skutek A. Ponieważ 
P(B; 0 A)= P(A|B;) : P(B;), więc otrzymujemy tzw. wzór Bayesa 


P(A|B.) P(B.) 


P(B;|A) = Z PAJB) P(B) 


(10) 


Wzór Bayesa nazywamy wzorem na prawdopodobieństwo przyczyny. 


e Przykład 2.4. 


Na stole leży 5 długopisów. Każde przypuszczenie (hipoteza) dotyczące liczby zepsutych 
długopisów jest jednakowo prawdopodobne. Wybrany losowo długopis okazał się zepsuty. 
Które przypuszczenie dotyczące liczby zepsutych długopisów jest najbardziej prawdopo- 
dobne? 

Rozwiązanie. 

Oznaczmy przez H; hipotezę, że liczba zepsutych długopisów wynosi 4, 4 = 0,1, 2,3,4,5. 
Mamy P(H,) = E. Niech A oznacza zdarzenie ” wylosowany długopis jest zepsuty”, wów- 
czas P(A|H,) = z. Zatem, ze wzoru na prawdopodobieństwo całkowite otrzymujemy 


a stąd 


poja U Bra 


Przy warunku, że wylosowano zepsuty długopis najbardziej prawdopodobne jest, że wszyst- 
kie długopisy są zepsute. 

W statystyce prawdopodobieństwa hipotez P(H;) nazywamy prawdopodobieństwami 

a priori (przed doświadczeniem), prawdopodobieństwa P(H,|A) - prawdopodobieństwa- 
mi a posteriori (po doświadczeniu). 


e Przykład 2.5. 
Telegraficzne przesyłanie informacji polega na wysyłaniu sygnałów: 0 albo 1. Przy przesy- 
łaniu 0 przekłamanie występuje w dwu przypadkach na trzydzieści, a przy przesyłaniu 1 
przekłamanie występuje w dwu przypadkach na czterdzieści. Stosunek liczby wysyłanych 
0 do liczby wysyłanych 1 wynosi 5 : 3. Obliczyć prawdopodobieństwo, że: 
a) wysłano 0, jeżeli wiadomo, że odebrano 0, 
b) wysłano 1, jeżeli wiadomo, że odebrano 1, 
c) wysłano 1, jeżeli wiadomo, że odebrano 0. 
Rozwiązanie. 
a) Oznaczmy przez Bı zdarzenie ” wysłano 0”, przez Bə zdarzenie ” wysłano 1”, przez 
A, zdarzenie odebrano 0” oraz przez A» zdarzenie odebrano 1” . Wiemy, że 


P(A»|Bi)=Ż, P(A;|Bo)=2Ż, 584 =5 


30 40” P(B2) 3: 


10 


Zgodnie ze wzorem Bayesa, prawdopodobieństwo zdarzenia ” wysłano 0, jeżeli odebrano 
0” równe jest 


P(B,|A,) = PAY PGW 


Zdarzenia Bı, Bə spełniają założenia twierdzenia o prawdopodobieństwie całkowitym, 
więc 
P(A) =P(A,|B;) : P(B) + P(4,|Bo) - P(B>). 


Ponieważ P(Ba) = $P(B4) oraz P(4,|By) = 1 — P(Ap|B;), więc 


= 1- P(A>|B1) — 280 ~ 
P(B;|44) = IPB) PANNE = żs9 © 0,969. 
b) Zgodnie ze wzorem Bayesa, prawdopodobieństwo zdarzenia ” wysłano 1, jeżeli odebrano 
1” równe jest 


A -P(B 
P(Bą| Aa) = Pa >) PB) c 2). 


Zdarzenia Bı, Bo spełniają założenia twierdzenia o prawdopodobieństwie całkowitym, 
więc 


P(A) = P(A>|B4) : P(B;) + P(A2|B2) : P(B2). 
Ponieważ P(B,) = 3P(B>) oraz P(A>|B>) = 1 — P(4,|B>), więc 
= 1-P(A1|B>) RE 
P(B2|A2) = PaB) PAB AR 


e Przykład 2.6. 
Przeciętnie 3% wyprodukowanych elementów ma wadę. Do wykrywania wady stosuje się 
test, który z prawdopodobieństwem 0,9 wskazuje wadę (wynik testu pozytywny), jeżeli 
element ma wadę i z prawdopodobieństwem 0, 95 nie wskazuje wady, jeżeli element jej nie 
ma. 
a) Jakie jest prawdopodobieństwo, że element ma wadę, jeżeli wynik testu jest pozytywny? 
b) Jakie jest powyższe prawdopodobieństwo, jeżeli element poddamy testowi dwukrotnie 
i za każdym razem otrzymamy pozytywny wynik testu? 
Rozwiązanie. 
a) Oznaczmy przez W zdarzenie ”element ma wadę” i przez N zdarzenie ” element nie ma 
wady.” Zdarzenia te są rozłączne, W UN = Q oraz 


P(W)=0.03, P(N) = 0.97. 


Niech D oznacza zdarzenie ” wynik testu jest pozytywny”. Z danych zawartych w zadaniu 
wynika, że 


P(D|W)=0.9, P(D'|N) = 0.95. 


Do obliczenia prawdopodobieństwa zdarzenia ”element ma wadę, jeżeli wynik testu był 
pozytywny”, czyli prawdopodobieństwa warunkowego P(W |D) wykorzystamy wzór Bay- 
esa. 

— PDIW):.P(W) _— P(DIW)-P(W) — 0.9-0.03 — 
P(W|D) = P(D) — P(DIW)-P(W)+P(DIN)-P(N) _ 0.9-0.03+0.05-0.97 0.358 


ponieważ z własności prawdopodobieństwa warunkowego wynika, że 
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P(DIN)=1-P(D'|N). 
Zatem, jeżeli wynik testu jest pozytywny, to prawdopodobieństwo, że losowo wybrany 
element ma wadę wynosi 35, 8%, czyli wśród elementów, które test wskazuje jako wadliwe 
tylko 35, 8% elementów ma wadę! Co wpływa na taką jakość testu? Jaki test byłby lepszy? 
Dla wad, które występują rzadko należałoby wykorzystywać testy o większych wartościach 
P(D|W) oraz P(D'|N). 
b) Obliczamy, jakie jest prawdopodobieństwo, że element ma wadę, jeżeli test przeprowa- 
dzony na nim dwukrotnie dał wyniki pozytywne. 
Niech A oznacza zdarzenie ” test przeprowadzony dwukrotnie na elemencie dał za każdym 
razem wynik pozytywny”. Wówczas 


P(AJW) = P(D|W).-P(DIW) oraz P(A|N) = P(D|N)- P(D|N). 


= P(A|W)-P(W) = (0.9)2.0.03 3 piba 


dwukrotnie zastosowany test dał wyniki pozytywne, to prawdopodobieństwo, że element 
ma wadę wynosi 0.909. Prawdopodobieństwo trafnej diagnozy znacznie wzrosło, gdy test 
przeprowadzilśmy dwukrotnie! tatwo policzyć, że prawdopodobieństwo zdarzenia ” element 
ma wadę, jeżeli test przeprowadzony trzykrotnie na tym elemencie dał za każdym razem 
wynik pozytywny” równe jest 0.994. 


e Przykład 2.7. 
Prawdopodobieństwo, że pogoda w danej miejscowości jest taka sama, jak dnia poprzed- 
niego równe jest a dla dnia deszczowego i b — dla dnia bezdeszczowego (a,be(0,1)). 
Prawdopodobieństwo, że pierwszy dzień roku jest deszczowy równe jest p4. Obliczyć praw- 
dopodobieństwo p,, że n-ty dzień roku jest deszczowy. 


Rozwiązanie. 


Oznaczmy przez D, zdarzenie "n-ty dzień roku jest deszczowy”, a przez B, zdarzenie 
”n-ty dzień roku jest bezdeszczowy”. Wówczas dla dowolnego n > 1 mamy: 


PDA DOZ: ZB BSK PDŁAJBJET=W PB Do i=u 

Policzmy 
po = P(D>) = P(D2|D1) - pr + P(D2|B;) ; (1 — pı) = api + (1 — b)(1 — pı) 

a+b—1)+1—b, 
D) = P(D3|D>) - pa + P(D3|Bo) - (1 — p2) = apa + (1 — b)(1 — p2) 
a+b—1)*+ (1 -—b)(a+b—1)+1—b, 
D4) = P(D4|D3) : pz + P(D4|E3) : (1 — ps) = pla +b- 1) +1—b 
= p(a+b-1)+ (1 —b)(a+b-— 1) + (1 —b)(a+b—1)+1—b, 


p = P 


pą = P 


Posługując się zasadą indukcji matematycznej można pokazać, że 


Pn = pila +b- 1) 1+ (1-b) Y (a+b-1)" 


k=0 
1—(a+b-— 1)" 
-= N E, ub waiza AS 
platba (1 
1— b) + (1—b)(a+b-— 1)! 
= R b-1) -————-— 
iec icz OPER G= a) + (1-5) 
Ponieważ a+b—1e(0,1), więc lim (a +b-— 1)" =0istąd lim M = oC: 
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2.3 Niezależność zdarzeń. 
Definicja 2.3. Dwa zdarzenia A i B nazywamy niezależnymi, jeżeli 
P(ANB) =P(A) -P(B). 


Jeżeli P(B) > 0, to z niezależności zdarzeń A i B wynika, że P(A|B) = P(A), czyli, jak 
się potocznie mówi, zajście zdarzenia B nie ma wpływu na prawdopodobieństwo zajścia 
zdarzenia A. Oczywiście dla dowolnego A € F zdarzenia A i Q są niezależne. Podobnie 
- zdarzenia A i Ø są niezależne. Jeżeli zdarzenia A i B są rozłączne i mają niezerowe 
prawdopodobieństwa, to nie mogą być niezależne. 


Zdarzenia Ay, Ao, A3,... nazywamy rodziną zdarzeń niezależnych, jeżeli dla każdej 
skończonej ilości zdarzeń A;,, A;,,..., A,, zachodzi równość 

P(A OQ- 0 Arn) =P(4,)::-P(4;,). (11) 
Fakt 2.4. Jeżeli zdarzenia A, Á2,..., An są niezależne, to niezależne są także zdarzenia 


Bı, By,..., Ba, gdzie B; = A; lub B; = A, dla i = 1,2,...,n. 


Przykład 2.8. 
Wyrazić prawdopodobieństwo sumy n niezależnych zdarzeń A;, i = 1,2,...,n za pomocą 
prawdopodobieństw poszczególnych składników. 
Rozwiązanie. 
Ponieważ 
UA=0Nf 4 


więc korzystając z własności prawdopodobieństwa i z niezależności zdarzeń 4,, i = 
1,2,...,n możemy napisać 
P(U A) =1-P(Ń A7=1- fi PA) = 1- [1 A — PA). 
i=1 i=1 i=1 i=1 


czyli 


P(U A) =1-1I[0 -P(4). (12) 


W tym momencie warto sobie przypomnieć, jak się oblicza prawdopodobieństwo sumy 
zdarzeń, o których nie wiadomo, czy są parami rozłączne (wzór ??). Widać, że, gdy 
zdarzenia A; są niezależne, ten dość skomplikowany wzór można mocno uprościć. 


Przykład 2.9. 


Trzech kontrolerów jakości pracuje niezależnie. Pierwszy wykrywa 90% wad, drugi - 80% 
a trzeci - 60%. Jaki procent wad wykrywają łącznie” Jaki procent wad wykrywa trzeci 
kontroler a nie wykrywa pierwszy ani drugi? 


Rozwiązanie. 
Niech A; oznacza zdarzenie ” wadę wykrył i-ty kontroler”. Wówczas: 


P(Ai) = 5, P(Az) =,  P(43) = g 


10? T0? 


Wada zostanie wykryta, gdy zajdzie zdarzenie A = A4 U Á U A3. Wykorzystując wzór 
na prawdopodobieństwo zdarzenia przeciwnego oraz prawa de Morgana otrzymujemy 
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P(A U Az U 43) = 1 — P(A} N 44 N 43). 


Ponieważ z niezależności zdarzeń 44, A>, Az wynika niezależność zdarzeń A, A3, Ag, więc 
P(A; U Ao U A3) = 1 — P(A4,)P(A,)P(A;) = 1 — 0,1-0,2- 0,4 = 0,992. 
Zatem łącznie kontrolerzy wykrywają 99,2% wad. 


Zdarzenie ” spośród trzech kontrolerów wadę wykrył tylko trzeci kontroler” można zapisać 
jako zdarzenie A; N A, N A3. Ponieważ te trzy zdarzenia (tzn. A, A>, A3) są niezależne, 
więc 


P(A, 0 A, N 43) = 0,1 -0,2-0,6 = 0,012. 


e Przykład 2.10. 
Niezawodnością urządzenia nazywamy prawdopodobieństwo tego, że będzie ono praco- 
wać poprawnie przez czas nie mniejszy niż T. Obliczyć niezawodność urządzeń, których 
schematy przedstawiają poniższe rysunki. Liczby pı, po,... oznaczają niezawodności po- 
szczególnych, niezależnie pracujących elementów. 


a) pı p2 NE Pn 
p1 
pa 
b) Pa 
p5 
p2 
pı Pe 
P3 p4 
C ) p7 


Rozwiązanie. 

Niech A; oznacza zdarzenie, że i-ty element pracuje poprawnie co najmniej przez czas T. 
Wtedy p; = P(A;). Niech p oznacza niezawodność urządzenia. 

a) Urządzenie pracuje niezawodnie wtedy i tylko wtedy, gdy niezawodny jest każdy ele- 
ment, czyli, gdy zajdzie zdarzenie A:N A>]. ..N An. Wykorzystując niezawodność zdarzeń 
A; otrzymujemy 

b) Urządzenie pracuje niezawodnie wtedy i tylko wtedy, gdy przynajmniej jeden element 
jest niezawodny, czyli, gdy zajdzie zdarzenie Aq U A2 U... U An. Stąd 


i wykorzystując prawa de Morgana oraz niezależność tych zdarzeń, otrzymujemy 


n 


1): P(A) -...- P(A„) =1— [] (1— ps) 


i=l 


p=1— P(A,04,0...0A,)=1—P(A 
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c) Urządzenie pracuje niezawodnie wtedy i tylko wtedy, gdy zajdzie zdarzenie A; i zda- 
rzenie A> U(A3] A4) i zdarzenie Az U AgU Az. Wykorzystując wzór na sumę zdarzeń oraz 
ich niezależność otrzymujemy 

= pı(p2 + P3 ` Pa — P2 ` P3 ` pa) (Ps + Pe + Pr — Ps P6 — Pe Pr — ps ` Pr + Ds pepr): 


e Przykład 2.11. 


Rozważamy rodziny posiadające n dzieci. Niech A oznacza zdarzenie, że rodzina ma dzie- 
ci obu płci, a B - rodzina ma przynajmniej jedną dziewczynkę. Czy zdarzenia A i B są 
niezależne? 

Rozwiązanie. 

Przyjmując, że dzieci w rodzinie uporządkowane są np. według starszeństwa, oznaczmy 
przez Q zbiór ciągów n-elementowych o elementach 0 (dziewczynka) i 1 (chłopiec). Wów- 
czas 


Zatem 


i równość nigdy nie zachodzi. 


e Przykład 2.12. 


Wkładamy losowo n ponumerowanych kul do n ponumerowanych szuflad. Jakie jest praw- 
dopodobieństwo p,, że przynajmniej jedna kula trafi do szuflady o tym samym numerze? 
Obliczyć „im Pn- 

Rozwiązanie. 

Niech A; oznacza zdarzenie ”i-ta kula wpadła do i-tej szuflady, i = 1,2,...,n . Każda 
kula wpada niezależnie do każdej z szuflad z tym samym prawdopodobieństwem równym 
L. Stąd, wykorzystując wzór 12, otrzymujemy 


n— oo 


Zatem lim p, = 1-— lim BL =]-e". 


3 Zmienne losowe jednowymiarowe. 


3.1 Definicja oraz rozkład i dystrybuanta zmiennej losowej. 


Matematyczny opis doświadczenia losowego wymaga sprecyzowania przestrzeni probabili- 
stycznej, z którą mamy do czynienia. Wprowadzimy teraz pojęcia, które pozwolą uprościć 
i stworzyć jednolity opis doświadczenia losowego. Zacznijmy od przykładu. 


Badamy jakość czterech wyprodukowanych elementów. Każdy z nich może być dobry 
(oznaczmy to przez 1) lub wadliwy (oznaczmy to przez 0). Mamy wówczas 


Q = £(Wi, W2, W3, wą) W; E £0, 1,4 = 1,2,3,4}. 
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F jest zbiorem wszystkich podzbiorów zbioru Q. Faktycznie, najczęściej interesuje nas 
tylko liczba elementów dobrych wśród sprawdzanych, czyli wynikowi (w1, wą, wz, wą) przy- 
pisujemy liczbę (w1 + w2 + w3 + wą). Określamy zatem na Q funkcję X o wartościach w R 
wzorem 

X (w1, w2, W3, w4) = W + w + W3 + w4. 


Wartościami funkcji X jest jedna z liczb ze zbioru 40, 1, 2,3, 4). 
Zdarzenie {(w1, wą, wz, wą) : X (w1, wą, wz, wą) = k} oznaczamy krótko (X = k). 


© Definicja 3.1. Zmienną losową nazywamy każdą funkcję X : Q — R taką, że dla 
dowolnego a € R zbiór (X < a) = {w EQ: X(w) < a} jest zdarzeniem losowym, czyli 
fweQ: X(v) < a) eF dla dowolnego a e R. 


W dalszym ciągu zapisujemy krótko 
(WED r X(w)<a}=(X <a). 


Z własności rodziny F wynika, że zdarzeniami losowymi są też wszystkie zbiory postaci: 
(X <a), (X > a), (X 2a), (a< X <b), (a< X <b), (a< X <b). 


Ze zmienną losową wiąże się pojęcie jej dystrybuanty. 


© Definicja 3.2. Dystrybuantą zmiennej losowej X: Q — R nazywamy funkcję 
Fx: R — [0,1] określoną wzorem: 


© Twierdzenie 3.1. Funkcja F : R — R jest dystrybuantą pewnej zmiennej losowej wtedy 
i tylko wtedy, gdy : 


e F jest niemalejąca, 
e lim F(xz)=0, lim F(x)=1. 
e F jest lewostronnie ciągła, 


Z dwu pierwszych warunków wynika, że dla każdego x € R prawdziwa jest nierówność 


0< F(x)< 1. 
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e Przykład 3.1. 


Czy można dobrać stałe a,b tak, by funkcja F(x) była dystrybuantą pewnej zmiennej 


losowej? 
a 


ae gdy x< -—l1, 
F(x)= 4 0,5 gdy —1<z<l, 
b(2—2) gdy z>1. 
Rozwiązanie. 


Sprawdźmy, czy istnieją stałe a,b, dla których F(x) spełnia założenia twierdzenia 3.1. 
F(x) jest lewostronnie ciągła. 

„im ae” = 0, więc a > 0. 

„m b(2 — >) = 2b = 1, więc b = Z. Aby F(x) była niemalejąca musi zachodzić warunek 
GE: SA. > oraz b = z. Zatem dla każdej pary liczb (a,b), gdzie 0 < a < $, b = > funkcja 
F(x) jest dystrybuantą pewnej zmiennej losowej. Tylko dla pary liczb (5,3) F(x) jest 
ciągła na R. 


Jeżeli potrafimy dla każdego podzbioru borelowskiego B określić prawdopodobieństwo, 
z jakim X przyjmuje wartości w zbiorze B, to mówimy, że został określony rozkład 
zmiennej losowej X: 


Px(B) = P({w : X(w) € B}). 


W dalszym ciągu oznaczamy krótko 
P(X € B) = Px(B). 


Zauważmy, że rozkład zmiennej losowej spełnia aksjomaty prawdopodobieństwa. 


Rozkład zmiennej losowej jest jednoznacznie wyznaczony przez jej dystrybuantę, co jest 
treścią następującego faktu. 


e Fakt 3.3. Prawdziwe są następujące równości: 
1.P(X>a)=1— Fy(a), 
2. P(a<X<b)=Fy(b) — Fx(a), 
3. P(X=a)=Fyxf(a+) — Fx(a) 
(stąd, jeżeli Fx jest ciągła w punkcie a, to P(X=a)=0), 
4. P(X < a) = Fy (a*), 
5. Jeżeli X jest typu ciągłego, to P(X = a) =0 dla każdego a € R. 


Wyróżniamy dwa zasadnicze typy zmiennych losowych: zmienne losowe typu skokowego i 
zmienne losowe typu ciągłego. 


© Definicja 3.4. Mówimy, że zmienna losowa X jest typu skokowego lub X jest zmien- 
ną dyskretną, jeżeli X przyjmuje skończenie lub co najwyżej przeliczalnie wiele wartości 
£i, i€ I przy czym P(X = w) = pi > 0 oraz > pę= l. 
jEI 


Rozkład prawdopodobieństwa zmiennej dyskretnej nazywa się często funkcją prawdo- 
podobieństwa i zapisuje w postaci 


((1:,p;) tel) 
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Dystrybuanta Fy : R — [0,1] zmiennej dyskretnej ma postać 


Fy(x)= P(X < x) = Pi. 


{i:xi<x} 
Jest to funkcja schodkowa, lewostronnie ciągła o skokach o wartości p; w punktach x;, i €T. 


e Przykład 3.2. 


Zorganizowano następującą grę : gracz wyciąga z talii 52 kart dwie karty (bez zwracania). 
Jeżeli są to dwa asy - gracz otrzymuje 20 punktów; jeżeli dwie figury (król, dama, walet) - 
gracz otrzymuje 10 punktów; w każdym pozostałych przypadków gracz traci dwa punkty. 
Znaleźć rozkład zmiennej losowej X oznaczającej wygraną gracza. Wyznaczyć i narysować 
dystrybuantę X. 


3.2 Najważniejsze rozkłady dyskretne 
e Rozkład zerojedynkowy . 


Mówimy, że zmienna losowa X ma rozkład zerojedynkowy z parametrem p (rozkład 
B(1,p)), jeżeli X przyjmuje tylko dwie wartości oznaczane przez 1 i 0 (nazywane od- 
powiednio sukcesem i porażką) oraz 


P(X=1)=p, P(X=0)=q=l=p gdzie pe(0,1). 


Typowymi przykładami zmiennych o rozkładzie zerojedynkowym są zmienne losowe, które 
opisują jakość wyrobu (dobry, wadliwy), pracę urządzea dwustanowych czy wynik gry 
(wygrana, przegrana). 


e Rozkład dwumianowy - rozkład Bernoulli'ego 


Eksperyment ze zmienną losową o rozkładzie B(1, p) powtarzamy niezależnie n razy. Niech 
X oznacza liczbę sukcesów w n powtórzeniach. Wówczas 


PIX == WZ dla GEWL z 


Mówimy, że zmienna losowa zdefiniowana wyżej ma rozkład Bernoulli'ego z parametrami 
(n, p). Piszemy wówczas, że zmienna losowa ma rozkład B(n, p). Łatwo sprawdzić, wyko- 
rzystując wzór na | dwumian Newtona, że nieujemne wartości P(X = k) sumują się do 1: 


> PR = k) = > (pt *=P+9"=L 


Zauważmy, że jeżeli X, (i = 1,2,...,n) są niezależnymi zmiennymi losowymi, z których 
każda ma rozkład B(1, p), to zmienna losowa Y =X,+X2+...+X, opisuje łączną liczbę 
sukcesów w tych n próbach, czyli ma właśnie rozkład Bernoulli'ego z parametrami n, p. 
Wartość ko, którą zmienna losowa dyskretna przyjmuje z największym prawdopodobień- 
stwem, nazywamy najbardziej prawdopodobną wartością X. Jeżeli X ma rozkład 
B(n,p), to 

7 l (n+ 1)p lub (n+1)p—1 gdy (n+1)pelN 

[(n+1)p] gdy (n+ 1)pg N. 
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e Przykład 3.3. 


W pewnym biurze zainstalowano 10 drukarek. Każda z drukarek pracuje niezależnie śred- 
nio przez 12 minut w ciągu jednej godziny. 

a) Jakie jest prawdopodobieństwo, że w losowo wybranej chwili będzie włączonych 7 dru- 
karek? co najmniej 7 drukarek? 

b) Jaka jest najbardziej prawdopodobna liczba drukarek włączonych w danej chwili? 
Rozwiązanie. 

Jeżeli drukarki pracują niezależnie średnio przez 12 minut w ciągu jednej godziny, to 
zmienna losowa X oznaczająca liczbę drukarek włączonych w danym momencie ma roz- 
kład Bernoulli'ego z parametrami n = 10, p = 2 = ł 

Zatem prawdopodobieństwo, że w losowo wybranej chwili będzie włączonych 7 druka- 


7 
rek wynosi P(X = 7) = (®) (4) (4 = 0.00079, a prawdopodobieństwo, że w losowo 
7745) 15 


wybranej chwili będzie włączonych co najmniej 7 drukarek równe jest 


8 (10) (1y* [41076 
P(X>7)= > CO G) E) =0,00086. 
Zatem, jeżeli zasilanie drukarek ustalone jest na poziomie dla sześciu drukarek, to prawdo- 
podobieństwo przeciążenia (czyli P(X > 7)) równe 0,00086. Czyli średnio przeciążenie ma 
miejsce w ciągu 86 minut na 100000 minut (1 minute na 1157 minut). Czy te rozważania 
mogą pomóc w ustaleniu poziomu zasilania? 


e Przykład 3.4. 
Prawdopodobieństwo prawidłowo wykonanej czynności przez pracownika wynosi 0.99. 
a) Jakie jest prawdopodobieństwo, że wszystkie spośród 100 takich samych, niezależnie 
wykonywanych, czynności zostaną wykonane prawidłowo? 
b) Jaka jest najbardziej prawdopodobna liczba czynności wykonanych prawidłowo? 
Rozważyć powyższe pytania, gdy czynność będzię powtarzana 199 razy. 
Rozwiązanie. 
a) Liczba prawidłowo wykonanych czynności wśród 100 niezależnych powtórzeń opisanej 
czynności jest zmienną losową o rozkładzie B(100, 0.99). Wszystkie czynności wykonane 
prawidłowo opisuje zdarzenie (X = 100), więc 


100 


P(X = 100) = ()(0.99)90(0.01)? zz 0.3660. 


b) Ponieważ dla zmiennej losowej X liczba (n + 1)p = 101 - 0.99 = 99.99 nie jest całko- 
wita, więc najbardziej prawdopodobną wartością X jest część całkowita z tej liczby, czyli 
99. Zatem zmienna losowa z największym prawdopodobieństwem przyjmuje wartość 99 i 
prawdopodobieństwo to wynosi 


P(X = 99) = (%4)(0.99) (0.01)! zz 0.3697. 


Niech zmienna losowa Y określa liczbę czynności wykonanych prawidłowo przy 199 po- 
wtórzeniach. Y ma rozkład B(199,0.99), więc 
P(X = 199) = Ea (0.99)! (0.01)? œ% 0.1353. 


199 


Dla zmiennej losowej Y liczba (n + 1)p = 200 : 0.99 = 198 jest całkowita, więc Y z naj- 
większym prawdopodobieństwem przyjmuje wartości 198 albo 197. Prawdopodobieństwo 
to wynosi 


P(X = 197) = P(Y = 198) = 0.2720. 


19 


Widać, że nawet jeżeli prawdopodobieństwo sukcesu w jednej próbie jest bardzo duże, 
to prawdopodobieństwo samych sukcesów oraz największa wartość prawdopodobieństwa 
maleje dość szybko wraz ze wzrostem liczby powtórzeń. 


e Przykład 3.5. 
Co jest bardziej prawdopodobne: wygrać z równorzędnym przeciwnikiem nie mniej niż 3 
partie z 4 partii, czy nie mniej niż 5 partii z 8 partii? 
Rozwiązanie. 
Zmienna losowa X określająca liczbę wygranych spotkań w czterech partiach ma rozkład 
B (4, 3), a zmienna losowa Y określająca liczbę wygranych spotkań w ośmiu partiach ma 
rozkład B (8, r). Otrzymujemy zatem 


Pa-3-0676'-4 70-05-0676" 


Zatem bardziej prawdopodobne jest wygranie dokładnie trzech spośród czterech partii. 
PXA>3=()G) G) +06) 6) =6 
Prs- OOOHOO OOTO OOOH 


czyli bardziej prawdopodobne jest wygranie przynajmniej pięciu spośród ośmiu partii niż 
przynajmniej trzech spośród czterech partii. 


e Rozkład Poissona z parametrem A 
Można udowodnić, że prawdziwe jest następujące twierdzenie. 


© Twierdzenie 3.2. (Poissona) 
Jeżeli (X,) jest ciągiem zmiennych losowych o rozkładzie dwumianowym B(n, pn), przy 
czym lim Npn = A, to dla każdego k € N zachodzi równość 


Mówimy, że zmienna losowa X ma rozkład Poissona z parametrem A, A > 0, jeżeli 


dą 


Oczywiście łatwo sprawdzić, że P(X = k) > O oraz 


dla k=0,1,2,.... 


OO 


P(X =k) = $ $ee SE A eh eA. 
=0 ; ! 


k k=0 k=0 


Przybliżanie rozkładu Bernoulli' ego rozkładem Poissona jest stosowane w przypadku, gdy 
n jest duże (n > 50) a p— małe tak, by np(1 — p) < 9. 


Najbardziej prawdopodobną wartością zmiennej lososwej o rozkładzie Poissona jest 


Joz AlubA—1 gdy AeN 
E A] gdy AgN. 
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Przykład 3.6. 

Po mieście jeździ 1000 samochodów. Prawdopodobieństwo, że samochód ulegnie uszko- 
dzeniu w ciągu jednej doby równe jest p = 0,002. Jakie jest prawdopodobieństwo, że co 
najmniej jeden samochód w ciągu doby ulegnie uszkodzeniu? Zakładamy, że samochody 
psują się niezależnie. Jaka jest najbardziej prawdopodobna liczba uszkodzonych samocho- 
dów? Ile miejsc należy przygotować na stacjach obsługi,by z prawdopodobieństwem 0, 95 
było wolne miejsce dla uszkodzonego samochodu? 

Rozwiązanie. 

Oczywiście możemy tu skorzystać ze schematu Bernoulli'ego, czyli 

P(X > 1)=1-P(X = 0) = 1 — ('47)(0.002)7(0.998)109 ~ 0.8649, jednak takie oblicze- 
nia są żmudne. 


Zgodnie z Twierdzeniem 3.2 (biorąc A = 1000 - 0,002 = 2) otrzymujemy 


P(X > 1)=1- P(X =0)=1— że”? z 0.8647. 


Ponieważ (n + 1)p = 2.002 g IN, więc najbardziej prawdopodobną liczbą uszkodzonych 
samochodów jest ko = 2. 

Aby odpowiedzieć na ostatnie pytanie należy znaleźć takie l, gdzie l jest liczbą miejsc w 
stacji obsługi, że P(X < l) > 0,95, czyli P(X>1)=P(X>1+1) <0,05. 
Posługując się tablicami rozkładu Poissona znajdujemy l+ 1 > 6, czyli l >5. 


Przykład 3.7. 
Wiadomo, że 1% produkowanych żarówek to braki. Obliczyć dokładnie i w przybliżeniu, 
prawdopodobieństwo, że: 
a) wśród losowo wybranych 100 żarówek nie ma ani jednej wybrakowanej, 
b) wśród losowo wybranych 100 żarówek są 2 wybrakowane, 
c) jaka jest minimalna liczba żarówek, które należy sprawdzić, by prawdopodobieństwo 
znalezienia złej żarówki było nie mniejsze niż 0,9. 
Rozwiązanie. 
Niech X oznacza liczbę wybrakowanych żarówek wśród 100 wylosowanych. 
0 100 0 
a) P(X =0)= (9) (75) (B) ~ei a~ 0,368 (A=100-7z=1 
b) P(X > 2)=1-— P(X =0)-— P(X = 1) x 0.2642 
c) Niech Y oznacza liczbę wybrakowanych żarówek wśród n wylosowanych. Y ma rozkład 
Bernoulli' ego B(n, 10): P(Y > 1)=1-— P(Y =0) 20,9, więc P(Y =0)< 0,1. Zatem 


szukamy n takiego, że P(Y = 0) = 3 (A (5) < 0.1. Stąd n > 230. 


Przykład 3.8. 


Śrubki są pakowane w pudełka po 100 sztuk. Prawdopodobieństwo, że śrubka jest wy- 
brakowana wynosi 0,01. Ile sztuk należałoby dodać do każdego pudełka, aby w celach 
marketingowych można było powiedzieć, że z prawdopodobieństwem nie mniejszym niż 
0,9 w każdym pudełku jest co najmniej 100 sztuk dobrych? 

Rozwiązanie. 

Niech X oznacza liczbę elementów wybrakowanych w pudełku. Zmienna losowa X ma 
rozkład Bernoulli'ego z parametrami n = 100, p = 0,01. Zatem 


P(X=0)=(%)pP( — p)! © = 0,366. 
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Jeżeli rozkład zmiennej X przybliżymy rozkładem Poisoona z parametrem 
A = 100 - 0,01 = 1, to 
P(X = 0) = e™! x 0,367. 
Dodanie do pudełka kilku (ko) elementów tylko nieznacznie zmieni parametr A (np. dla 


ko = 2 jest A = 1,02 oraz P(X =0)=e"* ~ 0,36059). Korzystając z przybliżonego 
rozkładu zmiennej X szukajmy więc najmniejszej liczby kg, dla której 


P(X < ko) ze *(1+4+$G+...+ gy) > 0,9. 


Przyjmując dla uproszczenia A = 1 i korzystając z tablic rozkładu Poissona otrzymujemy 
ko == 2, bo 


e7! (1+1+ 4) = 0,9196. 
Dokładniej — dla A = 1,02 mamy 


2! 


P(X < 2) = e71 (1 + 18 + 192) > 0,9295. 


W rzeczywistości prawdopodobieństwo znalezienia co najmniej 100 sztuk dobrych w pu- 
dełku zawierającym 102 śruby (czyli dla X o rozkładzie Bernoulliego z parametrami 
n = 102, p = 0,01) wynosi 


P(X < 2) = (0,99) + (0, 99)"®* - 0,01 - 1, 02 + #74% . (0, 99)" . (0,01)? = 0,9169. 


e Przykład 3.9. 


Liczba komputerów, które mogą być zarażone wirusem przez pewną sieć ma rozkład Po- 
issona z parametrem A. W każdym zarażonym komputerze wirus niezależnie uaktywnia 
się z prawdopodobieństwem p. Jakie jest prawdopodobieństwo, że wirus uaktywni się w 
m komputerach? 

Rozwiązanie. 

Niech zmienna losowa X oznacza liczbę zarażonych komputerów. X ma rozkład Poissona 
z parametrem A. Ponieważ zdarzenia (X = k), k=0,1,2,... są parami rozłączne oraz 


U (X=k)=Q9i P(X =k) = gon > 0, więc spełnione są założenia twierdzenia o 


prawdopodobieństwie całkowitym. Niech zmienna losowa Y oznacza liczbę komputerów, 
w których wirus uaktywni się. Mamy: 


P(Y =m) =P(Y=m)n U(Xx= 1)) - P(U (Y =m) n (X =4))). 
k=0 k=0 
Ponieważ zdarzenia (X = k), k = 0,1,2,... są parami rozłączne, więc zdarzenia (Y = 
m) N (X = k), k = 0,1,2,... są też parami rozłączne. Ponadto dla k < m jest (Y = 
m) N (X = k) =0. Zatem, wykorzystując twierdzenie o prawdopodobieństwie całkowitym 
otrzymujemy 


P(Y = m) = X P(Y = m) n (X =4)) = X P(Y = m|X =k) - P(X =4)). 


Zdarzenie (Y = m|X = k) dla k > m oznacza, że spośród k zarażonych komputerów 
wirus uaktywni się w m komputerach, w każdym z prawdopodobieństwem p (m sukcesów 
w k próbach), czyli 
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P(Y =m|X =k) = (;)p"(l- p)”. 


Podstawiając wartości prawdopodobieństw i dokonując elementarnych przekształceń, otrzy- 
mujemy 
OŚ k 
PY=m) =). (,)p"(1— p)" - e Xp 
=Ą m oo z l+m 
—e (>) 5u PD) 


(pa) 
10 z, 


e 9. 


e Rozkład geometryczny z parametrem p. 


Eksperyment ze zmienną losową o rozkładzie B(1, p) powtarzamy niezależnie dopóki po- 
jawi się sukces. Niech X oznacza numer próby, w której sukces pojawił się po raz pierwszy. 
Wówczas, kładąc q = 1 — p otrzymujemy 


P(X = k) =p", dla k=1,2... 


Mówimy, że zmienna losowa zdefiniowana wyżej ma rozkład geometryczny z parametrem 
p. Prawdopodobieństwa zdarzeń (X = k) są wyrazami ciągu geometrycznego i sumują się 
do 1: 


> P(X =k) = X p =p X i Ip hL 
k=0 k=1 k=l 


Zmienną o rozkładzie geometrycznym wygodnie jest interpretować jako czas oczekiwania 
na pierwszy sukces, tzn. liczbę powtórzeń eksperymentu poprzedzających eksperyment, 
w którym po raz pierwszy otrzymaliśmy sukces. 


© Fakt 3.5. Jeżeli zmienna losowa X ma rozkład geometryczny, to dla dowolnych liczb 
naturalnych ng, k zachodzi równość: 


P(X > no + klX > ne) = P(X > k). 


O zmiennej losowej spełniającej warunek z Faktu 3.5 mówimy, że ma tzw. własność 
braku pamięci. Ciekawszym jest fakt, że rozkład geometryczny jest jedynym rozkładem 
dyskretnym posiadającym własność braku pamięci. 


e Przykład 3.10. 
Prawdopodobieństwo, że danego dnia w miejscowości A latem świeci słońce jest stałe i 
równe p. Jakie jest prawdopodobieństwo, że jeszcze co najmniej przez 7 dni będzie piękna 
słoneczna pogoda, jeżeli już od dwu tygodni świeci słońce? 
Rozwiązanie. 
Przy założeniu, że pogoda w danym dniu nie zależy od pogody w dniach poprzednich, 
zmienna losowa X określająca liczbę kolejnych słonecznych dni ma rozkład geometryczny 
z parametrem p. Zatem 


P(X > 21|X > 14) = P(X > 20|X > 13) = P(X > 7) =(1—p).. 
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e Przykład 3.11. 
Główna wygrana w totolotku to prawidłowe skreślenie 6 liczb spśród 49. Jakie jest praw- 
dopodobieństwo głównej wygranej za 1001 razem, jeżeli przez 1000 razy nie było głównej 
wygranej? 
Rozwiązanie. 
Sześć liczb spośród 49-u można wybrać na E 3 sposobów (uporządkowanie liczb nie jest 


istotne). Wygranej sprzyja tylko jeden spośrod o ? jednakowo prawdopodobnych sposo- 
bów, czyli prawdopodobieństwo głównej wygranej p równe jest 


P= E sj = TRR , 


można powiedzieć, że jest bliskie jeden do czternastu milionów. 

Grę w ” Toto-lotka” powtarza się i w każdym powtórzeniu prawdopodobieństwo głównej 
wygranej równe jest wyżej obliczonemu p. Niech zmienne losowa X określa numer losowa- 
nia, w którym główna wygrana pojawi się po raz pierwszy. X ma rozkład geometryczny z 
parametrem p. Prawdopodobieństwo głównej wygranej za 1001 razem, jeżeli przez 1000 ra- 
zy nie było głównej wygranej to prawdopodobieństwo warunkowe P(X = 1001|X > 1000). 
Zatem 


P(X = 1001| X > 1000) a P(X=10010X>1000) z P(X=1001)- 


P(X>1000) P(X>1000) 
Ponieważ 
P(X > 1000)= 5, P(X=k)= ), q*"p=q". BE, = q"0, 
k=1001 k=1001 1 
więc 


1000 


P(X = 1001|X > 1000) = tm = p, 


czyli prawdopodobieństwo wygranej za 1001 razem, jeżeli nie wygraliśmy przez pierwsze 
1000 jest takie samo, jak prawdopodobieństwo wygranej za pierwszym razem. Na tym 
polega własność ” braku pamięci” rozkładu zmiennej losowej, a ma ją opisana zmienna 
losowa X. Jeżeli wiemy, że zmienna losowa przyjęła wartość większą niż n, to wszystkie 
następne wartości n-+k są przyjmowane z takimi samymi prawdopodobieństwami, z jakimi 
przyjmowane są wartości k. ” Przeszłość, jeżeli dotrwamy do chwili n, nie ma wpływu na 
” przyszłość”. 


e Rozkład Pascala z parametrami r,p. 


Eksperyment ze zmienną losową o rozkładzie B(1, p) powtarzamy niezależnie dopóki po- 
jawi się r sukcesów. Niech X oznacza numer próby, w której r-ty sukces pojawił się po 
raz pierwszy. Wówczas dla r > 1, q = 1 — p mamy 


k-1 
Beie e dlak=r,r+1,..., gdzie r>1,0<p<1. 


Mówimy, że zmienna losowa zdefiniowana wyżej ma rozkład Pascala z parametrami r, p. 
Zmienną o rozkładzie Pascala wygodnie jest interpretować jako czas oczekiwania na pierw- 
szy r-ty sukces. 
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e Przykład 3.12. 


Na ulicy stoi sprzedawca gazet. Każdy z mijających go przechodniów kupuje gazetę z 
prawdopodobieństwem p = E. Niech X oznacza ilość ludzi mijających go do momentu, 
gdy sprzeda 100 gazet. Znaleźć rozkład zmiennej X. 

Rozwiązanie. 


X ma rozkład Pascala z parametrami r = 100, p = 5, 


e Rozkład hipergeometryczny 


Z populacji składającej się z N elementów jednego rodzaju i M elementów drugiego 
rodzaju losujemy n elementów. Niech X oznacza liczbę elementów pierwszego rodzaju 
wśród wszystkich wylosowanych. Wówczas 


(w) (>) 
e lk n=k 
P(X =k) = o j 
Mówimy, że wyżej zdefiniowana zmienna losowa ma rozkład hipergeometryczny z para- 
metrami N, M,n, n< N,n < M. Definicja jest poprawna, bo można sprawdzić, że 


Ż Pix =k)= 5 — =1. 


k=0 n 


Przybliżenie rozkładem Poissona możemy stosować również w przypadku, gdy zmienna 
losowa ma rozkład hipergeometryczny, gdzie N + M jest duże a liczba < zz ` n mieści się 
w przedziale (0, 10). 


e Przykład 3.13. 


Zauważmy, że w Przykładzie 1.7 liczba ryb oznaczonych wśród 1000 wyłowionych (jeżeli 
w jeziorze jest N + 1000 ryb) jest zmienną losową o rozkładzie hipergeometrycznym z 
parametrami N, M = 1000, n = 1000. 


e Przykład 3.14. 
Spośród liczb 1,2,...,35 losujemy pięć liczb. Jakie jest prawdopodobienstwo, że będą 
wśród nich cztery mniejsze od 21? Jaka jest najbardziej prawdopodobna ilość liczb wylo- 
sowanych mniejszych od 21? Porównać wynik ze średnią z takich liczb. 
Rozwiązanie. 
W tym zadaniu mamy do czynienia ze zmienną losową (określającą ilość liczb mniej- 
szych od 21 wśród wszystkich wylosowanych) o rozkładzie hipergeometrycznym, gdzie 
N = 20, M = 15,n = 5, więc: 


Warunek P(X = k) < P(X = k +1) równoważny jest nierówności 


(2) (54) < (u) (5_(erv) 


E ię ) , 
5 5 


której rozwiązanie daje k < 3. 
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Przykład 3.15. 

Pudełko kulek potrzebnych do zmontowania łożyska zawiera 10 sztuk o dodatniej odchył- 
ce od nominalnego wymiaru średnicy i 15 sztuk - o ujemnej odchyłce. Do zmontowania 
łożyska potrzeba 6 kulek, z których co najwyżej 3 mogą mieć dodatnią odchyłkę od nomi- 
nalnego wymiaru średnicy. Jakie jest prawdopodobieństwo, że monterowi, który wybiera 
6 kulek losowo, uda się zmontować łożysko? 

Rozwiązanie. 

Niech X oznacza liczbę kulek wśród 6 wybranych o ddatniej odchyłce. Monterowi uda 
się zmontować łożysko, X < 3. X ma rozkład hipergeometryczny z parametrami: N = 
10, M = 15,n = 6. Zatem P(X <)=1-P(X=4)-P(X=5)-P(X=6). 


e Rozkład wielomianowy 


Rozkład dwumianowy możemy uogólnić na przypadek n powtarzanych niezależnych eks- 
perymentów, z których każdy może mieć jeden z k (k > 2) wyników. Niech p; oznacza 
prawdopodobieństwo realizacji wyniku i-tego rodzaju w każdej próbie, p; e(0,1), i = 
1,2,...,k, pı +po+...+pk = 1, zaś X; niech oznacza liczbę wyników i-tego rodzaju w 
n powtórzeniach. Wówczas 


n! 


P(X rS ni X= no Ak = ni) S= l i (M Po" ee 
nę.::no::...* Th: 


gdzie n,, i=1,2,...,k są liczbami naturalnymi oraz m + na +...+nk=n. 


Przykład 3.16. 

Jakie jest prawdopodobieństwo, że w sześciocyfrowym kodzie wystąpią trzy zera, dwie 

piątki i jedna ósemka? 

Rozwiązanie. 

Zmienne losowe X;, gdzie i = 0,1,2,...,9 oznaczają odpowiednio liczbę zer, jedynek, ..., 

dziewiątek w sześciocyfrowym kodzie. p; jest prawdopodobieństwem wylosowania jednej 

z dziesięciu cyfr, czyli p; = 105 41=0,1,2,...,9. Zatem 

P(X = 3,X5 = 2, Xs = 1,Xı = X = X; = Xı = Xs = X7 = X = Xo = 0) 
(LL — 

31211! (5) (5) (5) = 0.00006. 

Przykład 3.17. 

Po wstępnej kontroli technicznej 10% wyrobów oceniono jako dobre, 5% - jako wadliwe, 


a 25% zdecydowano poddać dalszej kontroli. Jakie jest prawdopodobieństwo, że wśród 10 
wylosowanych wyrobów jest 7 dobrych, 2 — wadliwe i 1 należy poddać dalszej kontroli? 


Rozwiązanie. 

Niech zmienne losowe X1, X2, X3 określają odpowiednio liczbę wyrobów dobrych, wadli- 
wych i przeznaczonych do dalszej kontroli wśród 10 wylosowanych. Prawdopodobieństwa 
dla poszczególnej jakości wyrobów wynoszą odpowiednio: p; = 0.7, pə = 0.05, pz = 0.25. 
Zatem 


P(X = 1, Xą = 2, X = 1) = =% (0.7)7 - (0.05)? - (0.25)! = 0.0185. 


W rozkładzie wielomianowym zmienna losowa X; określa liczbę elementów i-tego rodzaju 
wśród n elementów, więc X; ma rozkład dwumianowy B(n,p;). Zatem najbardziej praw- 
dopodobna liczba elementów przeznaczonych do dalszej kontroli spośród 10 wylosowanych 
to 2. 
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3.3 Zmienne losowe typu ciągłego 


© Definicja 3.6. Mówimy, że zmienna losowa X jest typu ciągłego, jeżeli istnieje nieujemna 
funkcja całkowalna fx :R —> R taka, że: 


Fx(x) = J fx(t)dt dla każdego x€ R. 


Funkcję fx nazywamy gęstością prawdopodobieństwa. 


Wiemy z analizy, że funkcja Fx jest wówczas ciągła. Ponadto - jest ona różniczkowalna 
we wszystkich punktach ciągłości funkcji fx i w punktach tych zachodzi równożć 


F(z) = fx(2). 
Nietrudno wykazać, że prawdziwa jest następująca charakteryzacja. 
© Twierdzenie 3.3. Funkcja f jest gęstością pewnej zmiennej losowej wtedy i tylko wtedy, 
gdy : 
e f(x) >0 dla każdego z € R, 


+00 
s | f(x)dr=1. 


e Przykład 3.18. 
Czy istnieje stała c, dla której podana funkcja jest funkcją gęstości: 


0, gdyz <0, 
f(x) = c, gdy0<z<l 
a, gdyxr>1 


Wykorzystamy twierdzenie 3.3. Dla funkcji f(x) mamy c > 0 oraz 


0 1 œ c 
i odz + | cda + | — dr = 2c=l1 
—00 0 1 a 


zatem dla c = > funkcja f(x) jest gęstością pewnej zmiennej losowej Y. Jej dystrybuanta 
ma postać 
0, gdyz <0, 
F(x) = iz, gdy0<x<1, 
1— +, gdyt>l1. 


Dg? 


Obliczmy, dla zmiennej losowej Y prawdopodobieństwo zdarzeń: 


11 
P(-2< X <1)= F(1)- F(-2)= | zdz= >, 
0 
P(O < X <4) = F(4) — F(0) = fra 277 
= Sa Ehk oe g 
oo 1 1 
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e Przykład 3.19. 


Dzienne zużycie energii (100 kWh =1) pewnej firmy jest zmienną losową X o gęstości: 


1(3+2x—2x7) dla O4z<3 
—_]) 9 f 
LOSE dla x <0 lub z >3. 


Jakie jest prawdopodobieństwo zdarzeń: (X > 3), (1< X <2)? 
Rozwiązanie. 


0 dla z < 0, 
Fx(x) = „Ue zT + ga” — ga —g dla 0<zr<3, 
1 dla x >3 


P(X > 4)=1-— Fx(})= $, P(1< X <2) = Fy(2) — Fx(1) = #. 
3.4 Najważniejsze rozkłady ciągłe 

e Rozkład jednostajny na odcinku [a,b]: 
Zmienna losowa X ma rozkład jednostajny na odcinku [a,b], jeżeli jej gęstość jest postaci 


z; dla x € [a,b], 


feta) = | 0, dla z ER) [at] 


Wówczas 
0 gdy t<Xa, 
Fx(t)=4 FE gdy t€ (a,b), 
1 gdy t >b. 
YA U 
1 1 
2 2 
o ż > O| ż > 
Rys. 1.03. Rys. 1.03. 
Gęstość rozkładu jednostajnego Dystrybuanta rozkładu jednostajnego 
na przedziale [0,2] na przedziale [0,2] 


e Przykład 3.20. 
Z przystanku autobusy odjeżdżają co 10 minut. Zakładamy, że rozkład T czasu przybycia 
pasażera na przystanek jest zmienną losową o rozkładzie jednostajnym. Obliczyć praw- 
dopodobieństwo, że pasażer będzie czekał co najmniej 4 minuty, mniej niż 3 minuty. 
Rozwiązanie. 


0 gdy £<0 
+ dla z € [0,10] | 
fr(z) = R PAR. Fr(t)= 4 ī gdy t€ (0,10], 
i R oóG 16 
0 dla z €R \ [0,10] 1 gdy t>10. 
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Zatem 


3 
P(T =F =—, 
(T <3) r(3) 10 
oraz ; 
P(T > 4) =1- P(T <4)=1- Fr(4) = 5 


e Przykład 3.21. 
Automat produkuje kulki metalowe o średnicy X będącej zmienną losową o gęstości 


= f5 dla xe[0.4,0.6], 
fx(2) a 0 dla ze R\ [0.4,0.6] 


Za zgodne z normą uznaje się kulki o średnicy z przedziału [0.41, 0.59]. Obliczyć praw- 
dopodobieństwo, że wybrana losowo z produkcji kulka spełnia wymagania normy. Jaka 
jest najbardziej prawdopodobna liczba kulek spełniających wymagania normy wśród 1000 
kulek. 

Rozwiązanie. 

P(0.41 < X < 0.59) = g = 0.9. Liczba Y kulek spełniających wymagania nor- 
my wśród 1000 kulek ma rozkład B(1000, 0.9), więc najbardziej prawdopodobna wartość 
zmiennej Y wynosi [1001 - 0.9] = 900. 


e Rozkład wykładniczy z parametrem A 


Zmienna losowa X ma rozkład wykładniczy z parametrem A, A > 0, jeżeli jej gęstość ma 
postać 


0 dla z < 0, 
jx(2) = je $ dla z> 0. 


Rys. 1.04. 


Gęstość rozkładu wykładniczego dla A= 1 


Wówczas dystrybuanta jest postaci 


0 dla z < 0, 
Fela) == | 1—e™ò dla z>0 


Rozkład wykładniczy posiada ” własność braku pamięci” przez co rozumiemy, że dla 
dowolnych nieujemnych z, s zachodzi równość: 


P(X>xz+s|X >s)=P(X > 2). 


Można także wykazać, że rozkład wykładniczy jest jedynym rozkładem ciągłym 
posiadającym własność ” braku pamięci”. 
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e Przykład 3.22. 


Czas pracy pewnego urządzenia jest zmienną losową X o rozkładzie wykładniczym z 
parametrem a = 107%. Wiadomo, że urządzenie pracowało 1000h. Jakie jest prawdopo- 
dobieństwo, że popracuje co najmniej 6000h? 

Rozwiązanie. 

Wykorzystując własność braku pamięci otrzymujemy 


P(X > 6000|X > 1000) = P(X > 5000+1000|X > 1000) = P(X > 5000) = 1—Fy(5000) 


e Rozkład normalny z parametrami m, o. 


Zmienna losowa X ma rozkład normalny z parametrami m, o (m, €R, o > 0), jeżeli 


jej gęstość ma postać 
2 
1 _ (c1-—m) l 


Wiemy już, że f fx(x)dz = 1. 
Rozkład ten oznaczać będziemy symbolem N (m, ø). 


Rozkład normalny jest najważniejszym rozkładem w teorii prawdopodobieństwa. Został 
wprowadzony w XVIIIw. przez Gaussa i Laplace’a Rozkład normalny, co niedługo przed- 
stawimy, stanowi dobre przybliżenie sum niezależnych zmiennych losowych. Z tego wzglę- 
du jest wykorzystywany do opisu losowych błędów pomiarów. Jeżeli błąd pomiaru nie- 
znanej wielkości jest sumą wielu małych losowych błędów dodatnich i ujemnych, to suma 
tych błędów ma rozkład bliski rozkładowi normalnemu. 


Dystrybuanty rozkładu normalnego N(m,o), czyli funkcji 


1 T _(t-m)? 
Fx(x) = a E 2? dt. 


nie można wyrazić przez funkcje elementarne. Wartości dystrybuanty rozkładu N(0, 1), 
czyli funkcji 


1 z t2 
p(s) = — j e`? dt. 
l ) V 21 J—oo 
podane są w tablicach. 


Wykres gęstości rozkładu N(0, 1) jest następujący 
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—0.5 


Rys. 1.04. 


Gęstość rozkładu N(0,1) 


Z symetrii wykresu gęstości względem osi Oy otrzymujemy wygodną w obliczeniach rów- 
NOŚĆ: 
Py(—z) =]- y(x). 


Okazuje się, że wartości dystrybuanty dowolnego rozkładu N(m,o) można obliczyć, 
znając wartości funkcji 9(z). 


p SE 
© Fakt 3.7. Jeżeli X ma rozkład N(m,o), to zmienna losowa Y = 2T ma rozkład 
o 


N(0,1) oraz 
Fx(1) = © (z -") 


Ostatni Fakt daje następujący, często wykorzystywany wzór 


Pa<x<)=0(+25) -o (12). 


Korzystając ze standaryzacji i z tablic rozkładu N (0,1) łatwo sprawdzić, że gdy X ma 
rozkład N(m,o), to 


P(m — 30 < X < m +30) = P(—3 < SĘ" < 3) = 26(3) — 1 = 2-0.9987 — 1 > 0.997. 


Oznacza to, że wartości zmiennej X z prawdopodobieństwem bliskim 1 zawarte są w 
przedziale (m — 30, m + 30). Własność tę nazywamy ” prawem trzech sigm”. 


e Przykład 3.23. 
Średnica metalowych kulek produkowanych przez automat jest zmienną losową X o roz- 
kładzie N(0.5,0.04). Za zgodne z normą uznaje się kulki o średnicy z przedziału [0.41, 
0.59]. Obliczyć prawdopodobieństwo, że wybrana losowo z produkcji kulka spełnia wyma- 
gania normy. Jaka jest najbardziej prawdopodobna liczba kulek spełniajacych wymagania 
normy wśród 1000 kulek? 
Rozwiązanie 
Szukane prawdopodobieństwo obliczymy dokonując standaryzacji zmiennej losowej X oraz 
wykorzystamy tablice rozkładu N(0, 1) 


0.41—05 X—05 _0.59—0.5 
0.04 0.04 0.04 
$(2.25) — ©(—2.25) = 20(2.25) — 1 = 0.9756 


P(0.41 < X < 0.59) = P( 


|= 


Zmienna losowa Y określająca liczbę kulek spełniających wymagania normy, wśród 1000 
kulek wyprodukowanych ma rozkład B(1000, 0.9756) i jej wartością najbardziej prawdo- 
podobną jest [1001 - 0,9756] = 976. 

Zauważmy, że obliczane w przykładzie wielkości są większe niż analogiczne wielkości dla 
zmiennej losowej opisanej w Przykładzie 3.20. 
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e Przykład 3.24. 
Czas sprawnej pracy drukarek pewnego typu (w dniach) ma rozkład N(1000,80).Jaki 
powinien być okres gwarancji, aby na 90% drukarka działała przynajmniej przez okres 
gwarancji ? 
Rozwiązanie 
Oznaczmy przez X czas sprawnej pracy drukarki, zaś przez T nieznany okres gwaran- 


cji.Zdarzenie (X > T) ma zachodzić z prawdopodobieństwem co najmniej 0.9, czyli 
P(X >T) > 0.9, ale 


X — 1000 T — 1000 T — 1000 T — 1000 
ROSNER) «jn. e pg Z, 
ORA p OBU > By 2221 0 © 

Szukamy T, dla którego 

T — 1000 
Es )>59 
E 


Z tablic rozkładu N(0, 1) odczytujemy, że -72% > 1.29. Ponieważ funkcja $(z) jest ro- 


snąca, więc T < 897.6. Zatem szukany okres gwarancji powinien mieć co najwyżej 897 dni. 
e Rozkład gamma z parametrami a,b > 0) 
Zmienna losowa X ma rozkład gamma z parametrami a,b > 0, jeżeli jej gęstość ma postać 
0 dla æ < 0, 


Ya „(z) == U a—1 —ba 
j l A 
NOW e dla z>0 


gdzie T (a) (tzw. funkcja gamma) zdefiniowana jest dla a > 0 za pomocą całki niewła- 
ściwej 


Tr(a) = Can 
0 


Całkowanie przez części daje wzór rekurencyjny 
T(a) =(a—1)7T(a—1) dla a>l, 
więc dla a= n € N otrzymujemy 
T(n) =(n- 1)! 


Dla naturalnych a = n rozkład gamma jest rozkładem sumy n niezależnych zmiennych 
losowych o rozkładzie wykładniczym z parametrem b. 


e Rozkład Cauchy'ego z parametrem A 


Zmienna losowa X ma rozkład Cauchy'ego z parametrem A, A > 0, jeżeli jej gęstość jest 
postaci 


A 
fx(z) = TETA 


Wówczas 


Fx(x) = 1 (arctg $ + =): 
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3.5 Funkcje zmiennych losowych 


Jeżeli g jest w miarę porządną (np. posiadającą co najwyżej przeliczalnie wiele punktów 
nieciągłości), funkcją zmiennej rzeczywistej, to dla dowolnej zmiennej losowej X: Q — R 
funkcja go X = g(X) : Q — R też jest zmienną losową. Mamy bowiem: 

(go X)7*((—00, a)) = X"*(g7'(-00,a)) € F. 

Dla zmiennej losowej X w pewnych przypadkach rozkład zmiennej losowej g(X), gdzie g 
jest funkcją zmiennej rzeczywistej można łatwo wyznaczyć. Jeżeli g jest Ściśle rosnąca, to 
dystrybuantę zmiennej Y = go X można wyznaczyć za pomocą dystrybuanty zmiennej 
X w następujący sposób: 


Fy(2) = P(Y < z) = P(g(X) < 2) = P(X <g (e) = Frig e). 


e Przykład 3.25. 
Miesięczny koszt prowadzenia przyzakładowego laboratorium jest zależny od liczby «x 
zatrudnionych w nim pracowników zgodnie ze wzorem y = 15000x + 100004/x. Koszty 
te traktujemy jako zmienną losową. Wyznaczyć funkcję prawdopodobieństwa zmiennej 
Y = 15000X + 10000v X, przyjmując następujący rozkład zmiennej losowej X: 


Di 


Rozwiązanie. 


Punkty skokowe y; zmiennej Y są postaci y; = 15000x; + 10000,/z;, co dla z; € £2,3,4,5) 
daje y; € 444142, 62321, 8000, 97361). Ponieważ g jest różnowartościowa, 


więc P(Y =y,) = P(X = x) = p, czyli otrzymujemy 


aeee ana 
Pi 


e Przykład 3.26. 


W ustalonym punkcie płaskiej folii znajduje się źródło promieniowania radioaktywnego 
wysyłające promienie równomiernie we wszystkich kierunkach. W odległości 1 od folii 
znajduje się równoległy do niej ekran, na którym obserwuje się błyski spowodowane pro- 
mieniowaniem. Niech X będzie zmienną losową oznaczjącą współrzędną punktu obserwo- 
wanego na ekranie. Korzystając z założenia, że kąt $(t) jest wartością zmiennej losowej © 
o rozkładzie jednostajnym na (0, r) wyznaczyć jej gęstość i dystrybuantę zmiennej losowej 
X Naszkicować obydwie funkcje. 

Rozwiązanie. 


Spójrzmy na rysunek. 


33 


A Niech (u, 1) będzie punktem, w 
którym umieszczone jest źródło 
promieniowania. Promienie wysy- 
łane są równomiernie we wszyst- 
kich kierunkach oznacza, że kąt 
Qt = 5 + w, gdzie tga = =: 
O > jest wartością zmiennej losowej o 
rozkładzie jednostajnym na (0, 7]. 


2 


Stąd 
Fx(t) = P(X <t) =P(© < g,) = P(* < 5 + arctg(t — u)) = 1|5 + arctg(t — u)| 


= į + rarctg(t — u), więc fx(t) = Fy(t) = ESC 
Otrzymaliśmy zatem ”przesunięty” rozkład Cauchy'ego. Jeżeli źródło promieniowania 
znajduje się w punkcie (0, 1) to mamy zdefiniowany wyżej rozkład Cauchy ego. 


© Definicja 3.8. Zmienne losowe X1, X2,..., Xn nazywają się niezależne, jeżeli dla do- 
wolnych ty,to,...,tn E R zachodzi równość 


P(X: <t, Xa < tz,..., Xn < tn) = |] | P(X: < ti) 


i=1 


e Przykład 3.27. 


Czas oczekiwania na połączenie w centrali telefonicznej dla każdego abonenta jest zmien- 
ną losową X o rozkładzie wykładniczym z parametrem A = 0,2s. Z centrali korzysta 
jednocześnie i niezależnie 100 abonentów. Obliczyć prawdopodobieństwo, że najkrótszy z 
czasów oczekiwania jest większy niż 5s a najdłuższy - mniejszy niż 10s. 
Rozwiązanie. 

Niech X, będzie czasem oczekiwania na połączenie i-tego abonenta, 

X, Xo, ...X190 Są zmiennymi losowymi niezależnymi. 


X = max| Xy, X, ... X100}, Y= minf X, Xp, ... X100}, 


Ponieważ 

(X < £) > y 1<k<100 (Xk < z), 
oraz 

(Y >x) 4> y 1<k<100 (Xk > x), 
więc 

Fx(z) = P(X <a) = P(X; < £, Xo <T, 222, X100 < w) 
= P(X < z)P(X2 < 2) ZE P(X100 < z) = (FG) 

Oraz 


Fy(x)= PY < x)= 1 — P(Y > x)= 1 — (1 — F(x))®™ 
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Dla zmiennych X4 o jednakowym rozkładzie wykładniczym otrzymujemy 


0 dla z < 0, _JO dla z < 0, 
PASS (l-e-8)0 da z>0. 00 = 1— e70 dla >O. 


Zatem; P(Y >5)=e7%, P(X <10)=(1—e2)0, 


3.6 Parametry rozkładu zmiennej losowej 


W praktyce istnieje potrzeba opisania zmiennej losowej przez podanie pewnych charakte- 
ryzujących ją liczb, zwanych parametrami rozkładu zmiennej losowej. 


Omówimy teraz najczęściej wykorzystywane parametry rozkładu zmiennej losowej. 


3.6.1 Wartość oczekiwana zmiennej losowej. 


© Definicja 3.9. Wartością oczekiwaną zmiennej losowej dyskretnej o rozkładzie 
{(xi, pi): i € I} nazywamy liczbę EX określoną wzorem: 


EX = 5, tipi, 
iEl 
przy założeniu, że suma © |x;|p; jest skończona. 
iEI 


Wartością oczekiwaną zmiennej losowej ciągłej o gęstości f(x) nazywamy liczbę EX okre- 
śloną wzorem: 


EX = ae 


—00 


-+00 
przy założeniu, że całka [| |z|f(x)dz jest skończona. 


Zauważmy, że wartość oczekiwana zmiennej losowej X jest odpowiednikiem znanego z 
fizyki pojęcia środka masy. W przypadku zmiennej dyskretnej prawdopodobieństwa p; 
interpretujemy jako masy skupione w punktach «x;, a przyjęty układ jednostek jest taki, 
że masa całkowita równa jest 1. W przypadku zmiennej ciągłej f(x) jest gęstością masy. 


Wprost z definicji (z własności szeregów i całek niewłaściwych zbieżnych) wynikają na- 
stępujące własności wartości oczekiwanej: 
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© Fakt 3.10. (Własności wartości oczekiwanej) 


RA 


E(aX +c) =aE(X)+c, dla a,c € R; 

Jeżeli P(X = ¢) = 1, to EX = c. W szczególności E(EX) = EX ; 

|EX| < E(|X|); 

Jeżeli P(X > 0)=1, to EX >20; 

Jeżeli Y = g(X), to EY = Xg(zi)pı w przypadku zmiennej dyskretnej oraz 


+00 
EY = f g(zx)f(x)da w przypadku zmiennej z gęstością f(x), o ile powyższy szereg 


i całka są zbieżne bezwzględnie; 


E(X+Y)=EX+EY; 


. Jeżeli X,Y są niezależne, to E(XY)=EX:EY. 


3.6.2  Wariancja zmiennej losowej. 


© Definicja 3.11. Wariancją zmiennej losowej nazywamy liczbę zdefiniowaną wzorem 


VarX = E(X — EX) 


Pierwiastek z wariancji, czyli D(X) = vVarX nazywamy odchyleniem standardo- 
wym albo dyspersją zmiennej losowej X, a wariancję oznaczamy często także sym- 
bolem D?(X). Wariancja pozwala ocenić, jak bardzo wartości zmiennej losowej różnią się 
od wartości oczekiwanej. Zaliczamy ją do grupy tzw. parametrów rozproszenia. Inter- 
pretując rozkład prawdopodobiedstwa jako rozkład masy jednostkowej (podobnie, jak w 
przypadku wartości oczekiwanej) widzimy, że wariancja jest odpowiednikiem występują- 
cego w fizyce pojęcia momentu bezwładności względem środka masy. 


© Fakt 3.12. Prawdziwa jest następująca równość 


VarX = EX? — (EXP. 


© Fakt 3.13. (Własności wariancji) 


1. 


VarX 20, 

VarX = 0 wtedy i tylko wtedy, gdy P(X = c) = 1 dla pewnej stałej c. 
Var(aX) = a? VarX, 

Var(X + c) = VarX, 


Jeżeli zmienne losowe X,Y są niezależne, to Var(X + Y) = VarX + Vary. 


Oczywiste jest uogólnienie dwu ostatnich własności Faktu 3.10 i ostatniej własności Faktu 
3.13 na przypadek dowolnej skończonej ilości zmiennych losowych. 
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Tabela wartości oczekiwanych i wariancji 


rozkład X 
dwumianowy B(n, p) p 


VarX 
(1-p) 


>a llt 
= 
| 


i 


Poissona z parametrem A 


geometryczny z parametrem p a 


S 
4 
— 
| 
2 
3 


= 
ko 
A 
— 
| 
> 
"3 


Pascala z parametrami r,p 


hipergeometryczny p=N(N+M)" np(1 — p) 


jednostajny na [a,b] 


R 
an 
g 
Zam 
o~ 
| 
Q 
— 
NM] 


wykładniczy z parametrem A 


x |bol 
= | 

-|| © 

+ 

= 
ii 
NIE I 
N) 


normalny N(m,o) 


gamma z parametrami a, b ab" * 


Q 
Si 
M 


e Przykład 3.28. 
Obliczyć wartość oczekiwaną i wariancję dla zmiennych losowych o następujących rozkła- 
dach: 
a) P(X = 49) = P(X = 51) = 4; 
b) P(Y = —100) = E, P(Y = 100) = 3, 
c) P(Z = 100) = P(Z = 0) = į}, P(Z =50)=4. 


Rozwiązanie. 


Opisane zmienne losowe są typu dyskretnego, więc: 

EX=49. +51- = 50, EY =—100- +100- ł = 50, EZ=100.++50-4+0- += 50. 
Te trzy zmienne losowe mają taką samą wartość oczekiwaną. Zwróćmy uwagę, że wartość 
oczekiwana zmiennej losowej zależy nie tylko od wartości, jakie ta zmienna losowa przyj- 
muje ale też od prawdopodobieństwa, z jakimi te wartości są przyjmowane. 

Wariancje tych zmiennych są następujące: 

VarX = (49 — 50)? - į + (51 — 50) -į = 1, 

VarY = (—100 — 50)? - + + (100 — 50) - + = 7500, 

VarZ = (100 — 50)? - + + (50 — 50)? - 4 + (0 — 50)? - 4 = 1250, 


Znajomość wariancji poprawia charakteryzację zmiennej losowej. Jeżeli wartości zmiennej 
losowej są odległe od wartości oczekiwanej, to wariancja jest duża. Mała wartość wariancji 
mówi, że wartości zmiennej losowej są bliskie jej wartości oczekiwanej. Zauważmy, że 
dla zmiennych losowych przyjmujących skończenie wiele wartości zawsze istnieje wartość 
oczekiwana i wariancja. 


e Przykład 3.29. 


Dla zmiennej losowej X, dla której EX = 2, VarX = 3, wykorzystując własności wartości 
oczekiwanej i wariancji, możemy obliczyć: 

a) E(4X — 9) = 4EX -9=—1, 

b) Var(—2X + 7) = (—2)*VarX = 12 

c) D(2X — 1) = /Var(2X — 1) = V2 VarX = 2V3. 
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e Przykład 3.30. 
Zmienne losowe X,Y są niezależne, X ma rozkład N(—3,2), Y ma rozkład N(4,3). Wy- 
korzystując własności wartości oczekiwanej i wariancji, możemy obliczyć: 
a) E(X —4Y +1) = 2EX -4EY +1 = —21, 
b) Var(X — Y) = VarX + (—1)*VarY = 5 
c) E(XY) = EX - EY = —12. 


e Przykład 3.31. 
Załóżmy, że niezależne zmienne losowe X; mają rozkład B(1,p), i = 1,2,...,n. Niech 
Y = 5 X,. Obliczyć EY oraz VarY. Dla jakiego p, VarY jest największa? 
=1 


1 
Rozwiązanie. 


Korzystając z własności wartości oczekiwanej mamy: 
EY = E(X Xi) =n-p. 

i=1 
Dzięki niezależności zmiennych X; możemy policzyć: 


VarY = © VarX; =n- p- (1-— p). 
i=1 

Widać, że VarY przyjmuje największą wartość dla p = z. 

Przypomnijmy też, że zmienna losowa Y ma rozkład B(n,p). Zatem otrzymaliśmy od 

razu parametry rozkładu B(n, p). 


e Przykład 3.32. 
X-EX 


VarX 


Pokazać, że jeżeli zmienna losowa X ma wariancję, to dla zmiennej losowej Y = 


zachodzą równości: EY = 0, Vary = 1. 

Rozwiązanie. 

Wykorzystując własności wartości oczekiwanej i wariancji zawarte w faktach 3.10, 3.13 
możemy policzyć 


oraz 


VarY = (x) Vax =Bxj= (zz) .VarX = 1. 


Przekształcenie powyższe nazywamy standaryzacją zmiennej losowej. 


3.6.3 Momenty. 


Wartość oczekiwana i wariancja są szczególnymi przypadkami parametrów rozkładu zwa- 
nych momentami. 


Momentem zwykłym rzędu k (k = 1,2,...) zmiennej losowej X nazywamy liczbę 
EX*. Momentem absolutnym rzędu k (k = 1,2,...) zmiennej losowej X nazywa- 
my liczbę E|X|*. Momentem centralnym rzędu k (k = 1,2,...) zmiennej losowej X 
nazywamy liczbę E(X — EX). 


© Fakt 3.14. Jeżeli istnieje moment absolutny rzędu k, to istnieją momenty absolutny, 
zwykły i centralny rzędów l < k. 
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Momenty centralne rzędu 3 i 4 wykorzystuje się do badania asymetrii rozkładu i stop- 
nia jego koncentracji wokół wartości oczekiwanej. 


© Definicja 3.15. Współczynnikiem asymetrii (skośności) nazywamy liczbę œs określoną 
wzorem 


E(X-EX)3 
D(X)3 


X3 = 


Współczynnik asymetrii dla rozkładu normalnego wynosi 0, a dla rozkładu wykładni- 
czego równy jest 2. 


3.6.4  Kwantyle. 


© Definicja 3.16. Kwantylem rzędu p, gdzie p € (0,1), rozkładu zmiennej losowej X na- 
zywamy każdą liczbę z, spełniającą warunek 


Fx(z,) < p < Fy(u,). 


Oznacza to, że P(X < z,) < pi jednocześnie P(X > xp) < 1-— p. Nie wszystkie 
kwantyle rozkładu zmiennej dyskretnej są jednoznacznie wyznaczone, a kwantyle rozkładu 
zmiennej typu ciągłego są jedynymi liczbami spełniającymi warunek Fy(x,) = p. 


Kwantyl rzędu 1 nazywamy medianą i oznaczamy m,. Dla zmiennej dyskretnej mamy 


zatem i 


T;LMEe MME 


a dla zmiennej typu ciągłego mediana spełnia równość 
F(me) = 0,5. 


1 3 ; 
Kwantyl rzędu 4 nazywamy kwartylem dolnym, a kwantyl rzędu 7 — kwartylem gór- 
nym. Mamy zatem 


P(uj < X sa) >i 


3 
4 
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3.7 Nierówność Czebyszewa i twierdzenia graniczne 


W tym rozdziale poznamy nierówności często wykorzystywane do szacowania prawdopo- 
dobieństw zdarzeń i zajmiemy się badaniem zbieżności ciągów sum zmiennych losowych. 


3.8 Ważne nierówności. 


Twierdzenie 1. (Nierówność Markowa) Jeżeli P(X ż0)=1, oraz EX < œ, to 
dla dowolnego k>0 zachodzi nierówność 


lub równoważnie: px 
P(X <k)>1- zm 
Dowód. 


Dla zmiennej X typu ciągłego o gęstości f(x) mamy np. 
Ex=| dz + | dr > kP(X > k), 
rf | af(ode>kPIX> b 


co daje żądaną nierówność. Analogiczny dowód można przeprowadzić dla zmiennej losowej 
dyskretnej. 


e Przykład 3.33. 


Śrubki są pakowane w pudełka po 100 sztuk. Prawdopodobieństwo, że śrubka jest wy- 
brakowana wynosi 0,01. Ile sztuk należałoby dodać do każdego pudełka, aby w celach 
marketingowych można było powiedzieć, że z prawdopodobieństwem nie mniejszym niż 
0,9 w każdym pudełku jest co najmniej 100 sztuk dobrych? 

Rozwiązanie. 

W rozdziale 3. rozwiązaliśmy to zadanie stosując przybliżenie rozkładu Bernoulli'ego roz- 
kładem Poissona. Teraz rozwiążemy je wykorzystując nierówność Markowa . Dla ustalone- 
go n niech X będzie zmienną losową określającą liczbę wybrakowanych śrubek w pudełku 
zawierającym n śrub. X ma rozkład Bernoulli' ego z parametrami n, p= 100: Chcemy, by 
P(n-X > 100) > 0.9 lub równoważnie P(n= X < 100) < 0.1. Wykorzystując nierówność 
Markowa dla £ = n — 100 możemy napisać 

15 

100* 


P(n — X < 100) = P(X > n — 100) < 


n: 
n— 


Wystarczy zatem znaleźć takie n, by zachodziła nierówność 


1 
Na 2 . . 2 2 
-iø Ś 0.1 równoważna nierówności n > 111. 


Parametrem, za pomocą którego można charakteryzować rozrzut wartości zmiennej loso- 
wej jest wariancja albo odchylenie standardowe. 


Rola wariancji jako miary rozrzutu zawarta jest w następującym twierdzeniu. 
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Twierdzenie 2. (Nierówność Czebyszewa) Jeżeli VarX < © , to dla dowolnego 
e>0 zachodzi nierówność: 


VarX 
aee ee E >. 
E 
lub równoważnie n 
a e E e 
E 


Dowód. 
Niech X będzie zmienną losową o skończonej wariancji. Wykorzystując nierówność Mar- 
kowa dla zmiennej Y = (X — EX }?, otrzymujemy: 

1 1 _ VarX 


P(X - EX| > £) = P(Y > æ) < ae ZEX —- EXP = a 


e Przykład 3.34. 


Zobaczmy, jakie oszacowanie daje nierówność Czebyszewa dla następujących zmiennych 
losowych: 


a) Zmienna losowa X ma rozkład normalny N(m,o). 
Ponieważ EX = m i VarX = o”, więc z nierówności Czebyszewa dla e = 30 mamy 


P(|X — m| > 30) < $ =}, 
podczas, gdy z prawa trzech sigm wynika, że prawdopodobieństwo tego zdarzenia jest 
mniejsze niż 0,003. 
b) Dla zmiennej losowej o rozkładzie prawdopodobieństwa 

P(X = —2) = P(X = 2) = 0,125, P(X =0)=0,75 
mamy 
EX = —2 . 0,125 +2- 0,125=0 oraz VarX =4-0,125+4-0,125=1 
Nierówność Czebyszewa z € = 2 ma postać 
P(|X — 0| > 2) < 

a jednocześnie prawdopodobieństwo tego zdarzenia wynosi 


P(|X| > 2) = P(X = —2 v X=2)= 


1 
4) 

czyli dla tej zmiennej w nierówności Czebyszewa zachodzi ” równość”. Zatem oszacowania, 
jakie daje nierówność Czebyszewa, nie można ” polepszyć”. 


e Przykład 3.35. 
W celu oszacowania nieznanej wadliwości p poddaje się kontroli kolejne elementy partii 
towaru. Niech X; = 0, jeżeli i-ty element jest dobry, X; = 1 - jeżeli i-ty element jest 
wadliwy, ż = 1,2,...,n oraz niech Zn = L > X,. Ile elementów należy skontrolować, by 


1=1 
P(|Z, — p| < 0.01) > 0,9? 
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Rozwiązanie. 


Obliczmy: 
EZ,=E (i 3 x, zie WARE Var( 
4-1 


nierówności Czebyszewa otrzymujemy: 


1x, = .n-p(l- p) = POP. Zatem z 


n 
i=l 


P(IZ, — p| < 0.01) > 1 — 
n 
Wyrażenie p(1 — p) przyjmuje wartość największą dla p = ż, więc 


p(p— 1) 1 
n- (0.01)? E7 (0.01)? 


Nierówność P(|Z, — p| < 0.01) > 0.9 zachodzi dla n spełniających nierówność 
1- mony > 09, czyli dla n > 25000. 


Przykład 3.36. 


Wykonujemy 80 rzutów kostką. Znaleźć najkrótszy przedział, w jaki z prawdopodobień- 
stwem nie mniejszym niż 0,9 wpada ilość otrzymanych szóstek. 

Rozwiązanie. 

Niech X będzie zmienną losową oznaczającą liczbę szóstek w 80 rzutach kostką. X ma roz- 
kład B(80, E). Z własności rozkładu Bernoulliego największe prawdopodobieństwa mają 
wartości zmiennej położone symetrycznie względem EX. Mamy: 

EX=€=$%, VarX =80:g:g= 4. 

Wykorzystując nierówność Czebyszewa otrzymujemy: 

P(X- 3|<e)>1— ge > 0,9. Stąd e > 10/10 = 10,54. Zatem najkrótszy przedział o 
szukanej własności ma postać (2.79, 23.87). 


3.9 Prawa wielkich liczb. 


Korzystając z nierówności Czebyszewa można udowodnić niektóre twierdzenia ważne dla 

zastosowań rachunku prawdopodobieństwa. 

Niech X, X2,... będą niezależnymi zmiennymi losowymi o wartościach oczekiwanych 
n n 

m, i wariancjach oż . Rozważmy Y, = >. X,. Wówczas VarY, = X o? i widać, że 
i=1 i=1 


rośnie ona wraz ze wzrostem n. Dlatego rozważa się Y, mnożone przez pewne liczbowe 


współczynniki A, , które ograniczają wariancję Y, . Najczęściej rozważa się A, = 1. 


Wtedy Y, jest średnią arytmetyczną zmiennych X1, Xa,..., X. 
Będziemy rozważać dwa typy twierdzeń granicznych: 


e Prawa wielkich liczb to twierdzenia, w których graniczna zmienna losowa ma 
rozkład jednopunktowy 


e Centralne twierdzenia graniczne to twierdzenia, w których graniczna zmienna 
losowa ma rozkład normalny. 


Niech X1, Xą,... będzie ciągiem zmiennych losowych o skończonej wartości oczekiwanej. 
Mówimy, że dla ciągu (X,) zachodzi: 


42 


e Słabe Prawo Wielkich Liczb, jeżeli 


p n 
Ve>0 im P (fo: Rx a>» BX Ż 
o n k=1 n k=1 


4)=" 


Ze Słabego Prawa Wielkich Liczb wynika, że z prawdopodobieństwem bliskim 1, wartośći 


n 
zmiennej losowej Ł 5X, różnią się dowolnie mało od swojej wartości oczekiwanej 
k=1 


Kiedy zachodzi Słabe Prawo Wielkich Liczb? 


Twierdzenie 3. Jeżeli X, są niezależne i „im. -zVar > x, =0, 


to dla ciągu (X,)zachodzi Słabe Prawo Wielkich | Liczb. 


D o w ó d. Korzystając z nierówności Czebyszewa możemy napisać: 


1 n 1 n n 
N k= N k= k21 o = 
Zauważmy, że jeżeli zmienne losowe X1, X3, ..., X, są niezależne i mają jednakowy rozkład 


o wartości oczekiwanej m i skończonej wariancji, to średnie Ł 5 Xk z prawdopodobień- 
k=1 

stwem bliskim 1 przybliżają m. 

Twierdzenie to jest prawdziwe również przy słabszych założeniach. 


e Przykład 3.37. 


W serii rzutów kostką niech 


1 gdy wypadnie 6, 

Xk = j ; 

0 w przeciwnym razie. 

Wówczas X; ma rozkład B(1, E). Z Prawa Wielkich Liczb wynika, że z prawdopodobień- 


n 
stwem bliskim 1, wartości zmiennej losowej t > X% różni się dowolnie mało od E. 
k=1 


e Przykład 3.38. 


Czy dla ciągu niezależnych zmiennych losowych o następujących rozkładach prawdopo- 
dobieństwa 


P(X, = Vn) = P(X, = -vn) = z: P(X, =0)=1- 3} 
zachodzi Słabe Prawo Wielkich Liczb? 
Rozwiązanie. 


Zmienne losowe X1, X>,... Są niezależne, ale nie mają jednakowych rozkładów. Sprawdzi- 
my zatem, czy spełniony jest warunek występujący w ostanim twierdzeniu. Ponieważ 


EX, = Vk- zz +(-Vk) zą+0=0 oraz VarXy=k: sg tk 


więc $5 VarX, = X Vk <n: yn, astąd 0< a Var (Xp Xp) < Ta = U; 
k=1 k=1 
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Zatem, dzięki twierdzeniu o trzech ciągach, 


„m 4 VAT Dk Xp = 0. 


Niech X1, X2,... będzie ciągiem zmiennych losowych o skończonej wartości oczeki- 
wanej. Mówimy, że dla ciągu (X,) zachodzi: 


e Mocne Prawo Wielkich Liczb, jeżeli prawdziwa jest równość 


n= n 


k=1 k=1 


© Twierdzenie 3.4. (Twierdzenie Chinczyna) 
Dla zmiennych niezależnych o tym samym rozkładzie i skończonej wartości oczekiwanej 


m zachodzi ER m 
P( lim AtA tirkama m) =ï 
noo n 


Pokażemy zastosowanie Mocnego Prawa Wielkich Liczb do obliczania całek oznaczonych, 
gdy wyznaczenie funkcji pierwotnej jest trudne lub wręcz - niemożliwe w sposób anali- 
tyczny. 

(Metoda Monte Carlo obliczania całek oznaczonych) 
Niech g(x) będzie funkcją całkowalną na przedziale |a,b| i oznaczmy przez I nieznaną 


b 
wartość całki | g(x)dx. Wiemy, że jeżeli X jest zmienną losową o rozkładzie jednostajnym 


na przedziale [a,b|, to 


b 
1 1 
Eg(X) = | sa) GRE "I. 
J b—a b—a 
Jeżeli X1, X2, ... są niezależnymi zmiennymi losowymi o rozkładzie jednostajnym na prze- 
dziale |a, b], to g(X1),g(X2),... są niezależnymi zmiennymi losowymi o jednakowym roz- 


1 


kładzie z wartością oczekiwaną —— 


- I. Zatem, na mocy Twierdzenia Chinczyna, mamy 


LĘ I 
P| lim =y g(X) = Eg(X,) = = 
(i g(Xr) = Eg(X1) =) Í 


k=1 


czyli z prawdopodobieństwem równym 1, dla dużych n, średnia z wartości 
g(X1),g(X2),...,g(Xn) jest dowolnie bliska liczbie = - I. Wystarczy więc wygenerować 
wartości zmiennych losowych X1, X2,..., X, z rozkładu jednostajnego na [a,b] i przyjąć 


n b 
liczbę =* © g(X;,) za oszacowanie całki I = f g(x)dz. 
k=l a 


Metody oparte na symulacji zmiennych losowych nazywają się metodami Monte Carlo. 
Można je również wykorzystywać do obliczania całek wielokrotnych. 
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3.10 "Twierdzenia Graniczne 
Mówimy, że dla ciągu (X,) zmiennych losowych o skończonych wariancjach zachodzi 


e Centralne Twierdzenia Graniczne, jeżeli dla dowolnego t € IR 


<t -Tdr = G(t) 


- j 
—— = —— e 
noo n Ą (12M 


n 
Xk 


n 
))XA-E 
k=1 1 


k 


Oznacza to, iż dystrybuanty zmiennych losowych Y, = dążą w każdym 


> Varx, 
punkcie t € R do dystrybuanty rozkładu normalnego N(0,1). 


Zauważmy, że Y, jest standaryzacją sumy Sn = >. X. 
k=1 


© Twierdzenie 3.5. (Twierdzenie Lindeberga-Levy'ego) 
Niech X, Xą,... będzie ciągiem niezależnych zmiennych losowych o jednakowym rozkła- 
dzie, wartości oczekiwanej m i wariancji 0 < o? < w . Wówczas dla ciągu (Xn) 
zachodzi Centralne Twierdzenia Graniczne. Zatem dla dowolnego telR prawdziwa jest 
równość ; 
X X4—nm 
lim P | EE — — <t| = G(i), 


noo oyn 
gdzie ®(t) jest dystrybuantą zmiennej losowej o rozkładzie N(0, 1). 


W praktyce Twierdzenie Lindeberga-Levy'ego wykorzystuje się przyjmując, że dla dużych 
n zachodzi przybliżona równość 


p(a< ZĘ <2) = (b) — (a), 


Stosując twierdzenie Lindeberga-Levy'ego do zmiennych o rozkładzie zerojedynkowym 
otrzymujemy: 


© Twierdzenie 3.6. (Twierdzenie Moivre’a-Laplace’a) 
Niech X, X2,... będzie ciągiem niezależnych zmiennych losowych o jednakowym rozkła- 
dzie B(1,p). Wówczas dla każdego t e R zachodzi równość 


lim P (> < r) = pt), 


gdzie S„=X1+X2+...+ X. 
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Ważne są uogólnienia Twierdzenie Lindeberga-Levy'ego idące w kierunku osłabiania za- 
łożeń tego twierdzenia. Okazuje się, że można zastąpić, dość restrykcyjne, założenie o 
jednakowym rozkładzie zmiennych X, założeniami dotyczącymi momentów trzeciego rzę- 


du. 


e Twierdzenie 3.7. (Twierdzenie Lapunowa) 


Jeżeli X, X»,... jest ciągiem niezależnych zmiennych losowych o skończonych trzecich 
momentach centralnych spełniającym warunek 


A> E|X, — EX,|* 
k=1 


lim ————— =(, 
noo n 
VarXy, 
1 


ze 
U 


to dla ciągu (X,) zachodzi Centralne Twierdzenie Graniczne. 


Istnienie trzecich momentów pozwala też oszacować dokładność przybliżenia w Central- 
nym Twierdzeniu Granicznym. 


e Twierdzenie 3.8. (Twierdzenie Berry-Essćena) 


Jeżeli X, X»,... jest ciągiem niezależnych zmiennych losowych o tym samym rozkładzie 
i E|X|$ < oo, to prawdziwa jest nierówność 


> X-E(X Xa) ECH. 


ai | o 0 E E) 
28 S VarX, vn(y VarXı) 


Komentarz. 


W prawach wielkich liczb widzieliśmy, że czynnik 4, = L (co daje po prostu śred- 


nią arytmetyczną zmiennych X, X2,..., X,.) jest tak silnie wygaszający, iż powoduje 
koncentrację rozkładu granicznego w jednym punkcie! W centralnych twierdzeniach gra- 
nicznych rozpatruje się współczynniki zmierzające wolniej do zera, mianowicie 4, = mał 


Wówczas okazuje się, że graniczna zmienna losowa ma rozkład normalny N(0,1). 


e Przykład 3.39. 


Z partii towaru o wadliwości 3% pobrano próbę 500-elementową. Obliczyć prawdopodo- 
bieństwo, że liczba elementów wadliwych w próbie nie przekroczy 4%. 
Rozwiązanie. 


500 
500 2 X; — 500 - 0,03 5 
PO X; < 20)= P CE E 


< —— | =6(041) = 0.6591. 
2 v500 -0,03-0,97 ` V145.5 G 
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e Przykład 3.40. 


Samolot zabiera 80 pasażerów. Przyjmując, że waga pasażera jest zmienną losową o roz- 
kładzie N(80,10) obliczyć prawdopodobieństwo, że łączna waga pasażerów przekroczy 
9000kg. 

Rozwiązanie. 

Niech X; oznacza wagę i-tego pasażera. Wówczas mamy 


80 
80 X X, — 80-80 
= 9000 — 6400 
P(X X, > 9000) = P | = 


El >. | a1- 0(29.07) £ 0. 
=] 10780 405 l ) 


e Przykład 3.41. 
Rzucamy 60 razy kostką. Znaleźć najkrótszy przedział, w jaki z prawdopodobieństwem 
większym niż 0,9 wpada ilość otrzymanych 'szóstek”. 
Rozwiązanie. 
Wykorzystując nierówności Czebyszewa (zad.3.4) otrzymaliśmy przedział [2, 79; 23,87]. 
Wykorzystując twierdzenia Moivre'a-Laplace'a otrzymujemy 


czyli ©(e) > 0.9 ie > 1.64. Zatem najkrótszy przedział jest postaci [7, 86; 18, 80]. 


e Przykład 3.42. 


Prawdopodobieństwo spostrzeżenia sputnika z ziemi z określonego punktu obserwacyjnego 
jest równe p = m przy każdym locie nad punktem obserwacyjnym. Znaleźć liczbę lotów, 
jaką powinien wykonać sputnik, aby z prawdopodobieństwem nie mniejszym niż 0,9 liczba 
X, spostrzeżeń sputnika była równa co najmniej 10. 

Rozwiązanie. 

Niech 


RE 1 gdy sputnik został dostrzeżony w k-tym locie, 
k — \ 0 gdy sputnik nie został dostrzeżony w k-tym locie. 


n 
X, € B(1, 5), Mk = 5 o? = pq 5% 1005 Yn 2 X; określa ilość spostrzeżeń 


sputnika w n lotach. 


Yan- _ 10-n-— 1 fa č g 
Pm >10) = P (py > TE) -S e zdł > 0,9. 
5VN 5VN V2n J- 


Mamy zatem: 1 — (a) =0,9 gdzie a = RR Z tablic odczytujemy, że a = —1,28. 


Zatem (po rozwiązaniu nierówności) otrzymujemy n > 147. 
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e Przykład 3.43. 
Komputer dodaje 1200 liczb rzeczywistych i każdą zaokrągla do najbliższej liczby cał- 
kowitej. Błędy zaokrągleń są niezależne i mają rozkład jednostajny na |[—4, 4]. Obliczyć 
prawdopodobieństwo, że błąd w obliczeniu sumy przekroczy 10. 
Rozwiązanie. 
Wprowadzając oznaczenia: X; - błąd zaokrąglenia i-tej liczby otrzymujemy: 


1200 


1 1 
EX; =0, oz = 15 oraz Var >. Xi = 1200: 75 = 100. 


Zatem 
1200 


1200 i i 10 
P X; >10) = P | &— > —— | = 1 — (1) = 0.1587 
> | vV100  v100 (U) 


e Przykład 3.44. 

Czas pracy diody ma rozkład wykładniczy o wartości oczekiwanej m = 900h. Zgroma- 
dzono zapas 100 diod. Obliczyć prawdopodobieństwo, że wystarczy ich na 100000 godzin 
pracy, jeżeli każdą diodę włączamy natychmiast po uszkodzeniu poprzedniej. Ile diod 
trzeba mieć w zapasie, by z prawdopodobieństwem większym niż 0,99 wystarczyło ich na 
99000 godzin pracy urządzenia. 

Rozwiązanie. 

Oznaczmy czas pracy i-tej diody przez X,. Wówczas X; ma rozkład wykładniczy z para- 
metrem A= 900: 


100 100 
P B Xi 2 100000) =1— P (o <` X< 100000) = 
i=1 


i=l 


100 X; — 90000 4 „10000 
10 - 900 10 - 900 


1-P (-10 < = 1 — (6(1,11) — 6(—10)) = 0.1335. 
N 

Szukamy N takiego,że P „ X;2 990000) > 0,99 czyli 
i=1 


N 
<),X; < 99000) < 0,01 
1 


< p PZE EZR a ZARA, 
900VN 900V N 


Wykorzystując CTG otrzymujemy 


N X- N -900 -N. 
p "m 99000 m 0.01 


(227 N -900 
900YVN 


albo równoważnie % (m) > 0.99 ~ (2.33). Szukane N spełnia nierówność 


9N — 9 - 2.33V N — 990 > 0, czyli N > 106. 


| < 0.01. 
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4 Zmienne losowe wielowymiarowe. 


4.1 Definicja i przykłady. 


© Definicja 4.1. Wektorem losowym n-wymiarowym (Zmienną losową n-wymiarową) nazy- 
wamy wektor n-wymiarowy, którego składowymi są zmienne losowe X; dlai=1,2,...,n, 


X (w)=(X1(w), Xo(w),..., Xz(w)) 


© Definicja 4.2. Dystrybuantą n-wymiarowej zmiennej losowej X nazywamy funkcję 
Fy (ti, t2,... tn): R” — R określoną wzorem 


Fy (ti t2,- tn) = P(X; < ti, Xə < lą,..., Xn < tn) 


Zajmiemy się bliżej zmiennymi losowymi dwuwymiarowymi. 


Dwuwymiarową zmienną losową (X,Y) przyjmującą co najwyżej przeliczalnie wiele war- 
tości {(x;, yj) :i E€ I,j e J} nazywamy dwuwymiarową zmienną losową dyskretną. 
Rozkład prawdopodobieństwa takiej zmiennej można przedstawić w postaci 


(((cuvs),pyu)), gdzie py = P(X = t, Y =y,), dla iel, jeJ. 


Dla zbiorów I, J skończonych wygodnie przedstawia się rozkład prawdopodobieństwa w 
postaci tabeli 


(AW | un |v || Uno 
lm |pn |p |-| Pin 
| Ge | Du | Poe |-| Pon 


| m _ | Bra | Pma |---| Pmn | 


Dystrybuanta takiej zmiennej jest funkcją schodkową 


F(x,y) = P(X <zx,Y <y)= 5 Pij. 


i, j;£i <T,Yj <y 


e Przykład 4.1. 


Rzucamy 3 razy monetą. Niech zmienna losowa X oznacza liczbę wyrzuconych orłów a 
zmienna losowa Y numer rzutu, w którym orzeł pojawił się po raz pierwszy. Łączny roz- 
kład prawdopodobieństwa wektora losowego X,Y opisuje następujaca tabela. 


UMI 
2 [05 [oi] 0 | 
[3 [os] o o] 
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Mówimy, że zmienna losowa (X,Y) jest typu ciągłego, jeżeli istnieje nieujemna funkcja 
całkowalna f(x,y) taka, że dystrybuanta ma postać 


ra = [ E 


W punktach ciągłości (xo, yo) funkcji f(x,y) 


oF 
Izy o W) = f (xo, yo). 


Dla borelowskiego zbioru A C R” mamy 


P((X,Y) € A) = J f(a,y))dzdy. 
A 


Następujące twierdzenie charakteryzuje dystrybuantę zmiennej losowej dwuwymiarowej 


© Twierdzenie 4.1. Funkcja F(x,y) jest dystrybuantą pewnej zmiennej losowej (X, Y )wtedy 
i tylko wtedy, gdy : 


e F(x,y) jest niemalejąca ze względu na każdą że zmiennych, 

e F(x,y) jest lewostronnie ciągła, 

e dla każdego x i każdegoy lim F(a,y)=0, lim F(a,y)=0 
1——00 y——00 


oraż zi F(x,y)=1. 


e dla każdych 11 < £2, Yı < %2 
F(x2,y2) — F (£1, Y2) — F (x2, y1) + F (21,41) > 0. 


Wnioskiem z twierdzenia 4.1 jest następująca charakteryzacja funkcji gęstości. 
© Twierdzenie 4.2. Funkcja f(x,y) jest gęstością rozkładu prawdopodobieństwa pewnego 
wektora losowego wtedy i tylko wtedy, gdy : 
e f(x,y) > 0 dla każdego (x,y) € R”, 
+00 +00 


« | | f(e y)dedy=1. 


—00 —00 


Znając rozkład prawdopodobieństwa wektora (X,Y) możemy wyznaczyć rozkłady praw- 
dopodobieństwa zmiennych X,Y. Nazywamy je rozkładami brzegowymi. W przypad- 
ku zmiennej losowej dwuwymiarowej dyskretnej (X, Y ) są one określone wzorami: 


pi = P(X EM) =X py, oraz pj = P(Y = y;) = X pi 
7 i 


Dla zmiennej dwuwymiarowej ciągłej (X,Y) tzw. gęstości brzegowe są następujące: 
x0=| fednwy frw) = | fade 


Rozkład wektora losowego (mówimy czasem rozkład łączny) wyznacza jednoznacznie 
rozkłady brzegowe, ale nie na odwrót. Rozkłady brzegowe wyznaczają rozkład łączny, gdy 
składowe wektora losowego są zmiennymi niezależnymi. 
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© Twierdzenie 4.3. Zmienne losowe X,Y są niezależne wtedy i tylko wtedy, gdy 
FKxwv)(r,y) = Fx(z) : Fy(y). 
W przypadku zmiennych dyskretnych warunek ten równoważny jest warunkowi 
Pik = Pipp dla wszystkich i,k 
a dla zmiennych typu ciągłego — warunkowi 


fæ y(t y) = fxlr)fy(y) dla wszystkich x,y € R. 


Oczywista jest postać powyższego twierdzenia w przypadku dowolnej skończonej ilości 
zmiennych losowych X1, X2,..., Xn- 


e Przykład 4.2. 
Zmienna losowa X jest liczbą spalonych zasilaczy w pracowni w ciagu dnia, zmienna loso- 
wa Y jest liczbą przepięć w sieci energetycznej. Łączny rozkład wektora losowego (X,Y) 
opisuje tabela 


a) Obliczyć P((X,Y) e £(2,0), (2, 1))). 

b) Wyznaczyć rozkłady brzegowe zmiennej losowej X oraz Y. Ile wynosi P(X = 1), 
P(Y =0). Obliczyć EX, EY. 

c) Czy zmienne losowe X, Y są niezależne? 


e Przykład 4.3. 
Wektor losowy (X,Y) ma rozkład o gęstości 


czy dla OXr<l O<y<yi 
0 poza tym 


f(z, y) -| 


a) Wyznaczyć stałą c. 
b) Wyznaczyć rozkłady brzegowe. 
c) Czy zmienne losowe X, Y są niezależne? 


e Przykład 4.4. 


z2 
Gęstość wektora losowego (X,Y) dana jest wzorem f(x,y) = Sa LEE 


a) Czy zmienne losowe X, Y są niezależne? 
b) Obliczyć P(X > 1). 
c) Obliczyć P((X,Y) € A), gdzie A = /(x,y): a” ry” < 1}. 
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4.2 Parametry rozkładu wektorów losowych 


Gdy dany jest rozkład wektora losowego (X,Y) oraz h : R” — R jest funkcją całkowalną, 
to dla Z = h(X,Y) 


J f hl(a,y)f(x,y)dzdy dla wektora losowego typu ciągłego 
EZ = Eh(X,Y) = | -> -% 
( ) D Alti, Yk)Pi,k dla wektora losowego typu dyskretnego 
i,k 


Dla wektora losowego (X,Y) kowariancją zmiennych X, Y nazywamy wielkość 
Cov( X,Y) = E(X — EX\(Y — EY) = EXY —EXFY. 
Jeżeli VarX > 0, VarY > 0, to definiujemy ważny parametr zwany współczynnikiem 


korelacji. 
Cov( X,Y) 


080 PaVa) 


© Twierdzenie 4.4. (Własności współczynnika korelacji): 
1. |p(X,Y| <1 
2.Jeżeli zmienne losowe są niezależne, to p(X,Y)=0. 
3. p(aX +b,cY + d) = sgn(ac)p(X,Y). 
4. p(X,Y) = +1 wtedy i tylko wtedy, gdy istnieją stałe a, b takie, że P(Y = aX +b) = 1. 


Współczynnik korelacji jest miarą zależności liniowej zmiennych X i Y. W przypadku, gdy 
p = 0, zmienne losowe nazywamy nieskorelowanymi. Jeżeli p(X,Y) = 0, to zmienne 
losowe moga być zależne. Swiadczy o tym poniższy przykład. 

e Przykład 4.5. 


Załóżmy, że zmienna losowa X ma rozkład N(0,o) i niech Y = X*. Można łatwo spraw- 
dzić, że Cov(X,Y ) = 0, a zmienne X,Y są zależne. 


e Przykład 4.6. 


Gęstość wektora losowego (X,Y) dana jest wzorem 


-ły cosg dla F<r<rm, OXy<2 


0 poza tym 


f(z, y) =| 


a) Znaleźć rozkłady brzegowe 
b) Wyznaczyć kowariancję oraz współczynnik korelacji zmiennych X,Y. Czy X,Y są nie- 
zależne? 


e Przykład 4.7. 


Współczynnik korelacji zmiennych losowych X, Y wynosi 0.25. Jaki współczynnik korelacji 
mają zmienne losowe 4X — 3 oraz —2Y +4? 
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